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A  leaky  waveguide  antenna  is  investigated  through  a  combination  of  the¬ 
oretical  analysis  and  numerical  simulation.  We  developed  a  design  proce¬ 
dure  based  on  the  analysis  of  Goldstone  and  Oliner,  for  an  aperture  in  the 
narrow  wall  of  a  rectangular  waveguide.  We  can  phase  scan  the  antenna 
by  adjusting  the  propagation  constant  of  the  guiding  structure,  and  we  can 
frequency  scan  it  by  taking  advantage  of  a  frequency  dispersive  behavior. 
A  combination  of  frequency  and  phase  scanning  can  be  used  to  steer  the 
beam.  We  describe  how  the  aperture  illumination  function  is  synthesized 
for  a  constant  width  aperture  and  we  present  the  design  equations.  We  use 
a  numerical  simulation  of  the  leaky  waveguide  section  to  obtain  the  radia¬ 
tion  efficiency.  We  then  calculate  the  realized  gain  to  evaluate  the  frequency 
scan  range  and  radiation  pattern  characteristics.  The  results  demonstrate 
that  the  main  beam  position  can  be  scanned  in  the  range  of  10°  to  30°  from 
broadside  over  a  narrow  frequency  range  without  corrupting  the  radiation 
pattern. 
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1.  Introduction 


A  traveling  wave  antenna  (TWA)  is  one  for  which  the  fields  and  currents 
producing  the  radiation  pattern  may  be  represented  by  traveling  waves. 
The  traveling  wave  phase  velocity  determines  the  angle  of  peak  radiation 
so  that  the  TWA  can  be  phase  scanned  by  the  adjustment  of  the  propaga¬ 
tion  constant  of  the  guiding  structure.  We  can  also  frequency  scan  the  TWA 
by  taking  advantage  of  the  frequency  dispersive  behavior  of  the  guiding 
structure.  Source  distributions  can  be  synthesized  to  produce  a  prescribed 
radiation  pattern  (e.g.,  pencil  beam,  sectorial,  etc),  and  a  combination  of 
frequency  and  phase  scanning  can  be  used  to  steer  the  beam.  TWA  design 
then  relates  the  source  distribution  to  the  phase  velocity  and  tries  to  imple¬ 
ment  the  desired  distribution  in  a  practical  structure.  A  large  scan  range 
is  desired  over  a  limited  (e.g.,  10  percent)  source  frequency  bandwidth. 
The  physical  structure  considered  is  a  rectangular  waveguide,  excited  by 
a  transverse  electric  (TE)io-mode  (or  Hio-mode),  with  a  long  slot  in  the 
waveguide  narrow  wall  radiating  into  a  half-space  (i.e.,  an  infinite  flange 
approximation).  The  radiating  aperture  can  be  represented  by  an  equiva¬ 
lent  magnetic  current  quasi-line  source  [1,  p  47]. 

The  perturbed  waveguide  supports  an  H-type  hybrid  mode  constructed  to 
satisfy  the  required  field  variation  as  described  in  appendix  A.  The  aper¬ 
ture  admittance  is  obtained  by  the  variational  technique  for  a  parallel  plate 
waveguide  radiating  into  a  half-space  as  shown  in  appendix  B.  The  equiv¬ 
alent  network  model  includes  an  internal  susceptance  to  account  for  the 
stored  energy  attributable  to  the  aperture  perturbation  on  the  waveguide 
[2].  A  transmission  line  description  in  the  transverse  direction  is  then  termi¬ 
nated  in  this  equivalent  admittance.  We  use  the  transverse  resonance  tech¬ 
nique,  as  described  in  appendix  C,  to  determine  the  transverse  wavenum¬ 
ber  and  thus  the  perturbed  waveguide  complex  propagation  constant.  The 
slotted  waveguide  section  then  is  a  nonuniform,  lossy  transmission  line  in 
which  the  energy  is  radiated  along  the  transmission  line  direction.  The  fi¬ 
nite  element  method  (FEM)  is  used  to  solve  the  nonuniform  transmission 
line  equations  for  the  transmission  and  reflection  coefficient  of  the  slotted 
section  from  which  we  calculate  the  antenna  efficiency.  This  network  model 
characterizes  the  aperture  illumination  function  and  we  obtain  the  radia¬ 
tion  pattern  by  direct  integration  of  the  source  current.  The  basic  numerical 
routines  to  produce  these  calculations  are  included  as  appendix  D.  A  thor¬ 
ough  understanding  of  the  radiation  by  this  type  of  dispersive  waveguide . 
is  required  to  develop  a  design  procedure  for  electronically  steered  TWA 
arrays  fabricated  in  a  rectangular  waveguide. 
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2.  Leaky  Wave  Antenna 


We  consider  a  thin,  rectangular  aperture  of  width,  d,  and  length,  L,  in  a 
terminated  rectangular  waveguide  having  width,  a,  and  height,  b,  shown 
in  figure  la  as  a  leaky  wave  antenna  structure.  For  narrow  apertures,  the 
waveguide  perturbation  is  small  or  rather,  the  aperture  is  loosely  coupled, 
so  that  the  phase  per  unit  length  is  slowly  varying  with  position  along  the 
guiding  structure.  The  leaky  waveguide  continuously  loses  energy  along 
its  length  because  of  radiation,  as  shown  in  figure  lb.  At  a  particular  op¬ 
erating  frequency,  /o  =  c/  Ao,  the  main  beam  appears  at  an  elevation  an¬ 
gle,  sin  Oq  —  (Ao/A g),  from  broadside,  depending  on  the  perturbed  guide 
wavelength,  Xg  (see  fig.  la).  The  Hm-mode  power  available  from  the  in¬ 
put  waveguide,  Pmc,  is  partially  reflected  by  the  slotted  section  with  re¬ 
flection  coefficient,  T in-  Then  the  input  power  to  the  TWA  depends  on  the 
input  impedance  mismatch,  Pm  =  Pj110  -  Preflect  =  Pi»c(l  -  |rio|2)-  We  de¬ 
fine  the  input  efficiency,  ?/i„  =  (1  —  |T io|2)^  as  a  measure  of  the  power  lost 
because  of  reflection,  since  Pm  =  r?jnPinc,  is  the  power  fed  into  the  slot¬ 
ted  section  (see  fig.  lb).  The  radiated  power  depends  on  the  fraction  of 
power  transmitted,  |7m|2/  through  the  slotted  section  (i.e.,  absorbed  in  the 
matched  load).  If  we  assume  that  the  power  not  transmitted  is  radiated, 
then  Prac i  =  Pjn  —  Ptrans  =  Pin(l  —  |Pio|2)-  Therefore,  the  radiation  effi¬ 
ciency,  i]r  =  (1  —  |Pio|2),  depends  on  the  power  lost  in  the  load  impedance. 
We  group  other  loss  mechanisms  into  a  single  absorption  coefficient.  Am, 
and  define  an  efficiency  as  ??c  =  (1  -  |Am|2).  The  total  antenna  efficiency, 
V  —  riinflrVcr  is  a  measure  of  the  power  lost  because  of  input  reflection,  dis¬ 
sipation  in  the  load  impedance,  and  other  losses,  or  Praf j  =  7?Pjnc. 


For  a  given  available  power,  Pmc,  corresponding  to  the  input  power,  Pin, 
in  the  leaky  waveguide,  we  denote  by  P(y )  the  power  radiated  per  unit 
aperture  length.  By  conservation  of  energy,  the  total  power  flowing  at  some 
point,  Po(?y),  within  the  waveguide  extending  from  y0  =  —L/2  to  y  =  L/2 
can  be  written  as 


po(y)  =  Pm  -  [VP(m- 

Jyo 


(1) 


Differentiating  equation  (1)  with  respect  to  y  and  defining  the  attenuation 
per  unit  length  by 


My) 


P{y) 

2Po(y)  ’ 


(2) 


one  finds  that  Po(y)  must  satisfy  the  differential  equation 


+  2a, •(?/))  My)  =  0  .  (3) 

A  wave  propagating  down  the  guide  then  is  continuously  losing  energy 
through  radiation  [3].  We  obtain  the  solution  to  equation  (3)  using  the  initial 
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Figure  1.  (a)  Slotted 
waveguide  leaky  wave 
antenna  oriented  along 
y- axis  (not  to  scale)  and 
(b)  general  leaky  wave 
coupling  structure 
showing  mail-beam 
angle. 


condition,  P[n  =  Po(yo)-  The  result  is 

Po(y)  =  -Pin  exp  (-2  f  ar(£)d£j  .  (4) 

Assuming  no  other  losses  (i.e.,  Aw  —  0  or  i)c  —  1),  the  fraction  of  power 
not  radiated,  q,  will  be  completely  absorbed  by  the  termination  at  y  =  L/2. 
Thus,  we  define 


(5) 


The  radiated  power  per  unit  length  is  proportional  to  the  squared  magni¬ 
tude  of  the  aperture  illumination  function,  f2(y),  which  is  proportional  to 
the  aperture  electric  (E-)  field.  That  is,  P{y)  =  (  f 2(y),  in  which  Q  is  a  con¬ 
stant.  Using  equation  (5),  we  find  that  this  proportionality  constant  is 


^  _  (1  q)  Pin  _  VrPin 

~W2 


with  I  = 


(6) 


in  which  jjr  =  1  -  q  is  equal  to  the  previously  defined  radiation  efficiency. 
Once  the  magnitude  of  the  aperture  illumination  function,  |f(y)|,  is  chosen 
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Figure  2.  Normalized 
attenuation  (in  Nepers) 
for  a  uniform  aperture 
illumination  function. 


for  the  desired  radiation  pattern,  we  can  calculate  the  required  attenuation 
function  for  the  leaky  wave  antenna  by  substituting  £f 2(y)  into  equation  (2) 
and  using  equation  (6)  to  find 


or(y) 


Vr 


(7) 


This  is  an  explicit  solution  for  ar(y )  in  terms  of  the  aperture  illumination 
function,  which  is  chosen  to  obtain  the  desired  radiation  pattern  charac¬ 
teristics.  Then  equation  (7)  determines  the  attenuation  function,  which  is 
synthesized  by  means  of  an  appropriate  coupling  structure.  The  fraction 
of  power  absorbed  in  the  load  can  be  made  small  but  not  completely  zero 
since  this  would  require  an  infinitely  long  antenna  or  an  ar(y)  too  large  to 
be  realized  by  physical  structures. 


For  uniform  illumination  with  j/0  =  -L/2,  equation  (7)  takes  a  simple  form 


ttr  (y) 


L_ 

nr 


1 

2 


-  (y  -  yo) ' 


(8) 


We  normalize  the  aperture  dimensions  and  the  results  to  the  free-space 
wavelength,  A0.  Using  equation  (8),  the  normalized  attenuation,  |ar(y)Ao| 
(in  Nepers)  is  shown  in  figure  2  for  L  =  20Ao,  a  uniform  aperture  illumina¬ 
tion,  and  the  radiation  efficiency,  t]r  =  70,  80,  and  90-percent  (—1.5,  -1  and 
—0.5  dB),  as  a  parameter.  To  obtain  a  uniform  illumination,  the  correspond¬ 
ing  attenuation  per  unit  length  would  require  an  aperture  with  flared  width 
profile  so  that  the  aperture  width  would  increase  with  distance  toward  the 
load.  For  a  fixed  aperture  length,  a  larger  width  or  flare  per  unit  length  in¬ 
creases  the  radiation  efficiency  but  does  not  significantly  increase  the  input 
reflection  coefficient.  That  is,  the  total  antenna  efficiency  is  approximately 
equal  to  the  radiation  efficiency  shown  in  figure  2. 
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Conversely,  for  a  specified  attenuation  function,  we  can  determine  the  aper¬ 
ture  illumination  function  as  follows.  Define  X(y)  =  f 2(y),  and  rewrite 
equation  (7)  as 


(9a) 

which  can  be  written  as 

X(y)=2ar(y){^-  fl(^)  . 

1 Vr  jyo  J 

(9b) 

Differentiating  equation  (9b)  with  respect  to  y  and  using  equation  (9a)  gives 

dy  =x(v){l -^or(y)-^r(y))  . 

(10) 

Solving  equation  (10)  gives 

X{y)  =  cT{ylar{y)eXp{  2 ■ 

(ID 

Then  in  terms  of  ar(y),  the  magnitude  of  the  aperture  illumination  function 
is 

\f{y)\  =  \JCar(y)  exp  (  -  J  ar{^)d^j  (12) 

in  which  C  —  X(yo)/ar(yo)  is  an  arbitrary  constant. 

For  constant  aperture  width  (normalized  to  waveguide  height),  d  =  b/ 30, 
the  attenuation  is  calculated  as  developed  in  the  next  section.  The  corre¬ 
sponding  normalized  aperture  illumination  function,  f(y)/f(y0),  is  shown 
in  figure  3,  which  indicates  that  for  constant  aperture  width,  the  illumina¬ 
tion  function  decreases  along  the  aperture  length.  The  results  correspond 
to  a  waveguide  with  ka  =  4.25  in  which  k  =  2-k/Xq  is  the  wavenumber.  For 
a  narrow  aperture,  the  aperture  length  must  be  significantly  increased  to 
obtain  r/r  >  90  percent.  For  a  given  aperture  width  profile,  the  aperture  il¬ 
lumination  function  is  determined  and  the  variation  in  the  waveguide  com¬ 
plex  propagation  constant  is  calculated  with  an  equivalent  network  model. 
The  basic  leaky  wave  antenna  design  then  simplifies  to  implementing  the 
required  ar  (y)  in  a  slotted  waveguide  section.  These  leaky  waveguide  an¬ 
tennas  can  be  readily  combined  to  form  a  two-dimensional  planar  array  by 
accounting  for  mutual  interactions  [4]. 

We  determine  the  attenuation  function  for  a  section  of  leaky  waveguide 
from  the  aperture  geometry  with  perturbed  waveguide  propagation  con¬ 
stant,  7  —  ar +j/3.  In  general,  the  attenuation  constant  is  not  independent  of 
the  phase  constant,  but  the  waveguide  cross  section  can  be  adjusted  to  keep 
(3  -  2ir/Xg  approximately  constant  while  ar  [2,  p  190]  is  varied.  Radiating 
elements  fabricated  in  the  waveguide  narrow  wall  allow  the  use  of  spe¬ 
cialized  waveguide  types  such  as  dielectric  loaded  or  ridged  waveguide, 
which  provide  additional  design  flexibility.  Typical  examples  are  shown  in 
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Figure  3.  Aperture 
illumination  function 
for  a  constant  width 
aperture  with  radiation 
efficiency  as  a 
parameter. 


figure  4a  for  radiating  elements  in  the  narrow  wall  of  a  rectangular  wave¬ 
guide  [1,  p  282],  We  have  investigated  tapered  and  flared  transition  sec¬ 
tions  as  shown  in  figure  4b  as  a  "discontinuity  minimizer"  to  avoid  exciting 
higher  order  modes  [1,  p  284]. 

As  an  example,  we  consider  an  aperture  with  a  constant  width  section, 
do  =  6/40,  and  total  length  (normalized  to  waveguide  width),  L  =  30a.  We 
include  either  a  tapered  input  section  or  a  flared  output  section  of  length 
L\  =  L/ 2  and  L 2  =  L/2,  respectively.  For  narrow  apertures,  the  input  re¬ 
flection  is  small  with  r)in  >  97-percent  (-0.1  dB).  Nevertheless,  the  phase 
variation  over  the  tapered  section  produced  an  antenna  pattern  with  in¬ 
creased  side  lobe  level  (SLL).  Although  not  described  in  detail,  the  results 
for  variable  aperture  width  demonstrated  that  tapered  or  flared  sections 
(see  fig.  4b)  are  undesirable  in  terms  of  SLL.  A  flared  output  section  pro¬ 
vides  a  better  approximation  to  a  uniform  aperture  illumination,  as  can 
be  surmised  from  equation  (12).  However,  the  flared  (or  tapered)  aperture 
width  corrupts  the  radiation  pattern,  as  shown  in  figure  5.  The  tapered  in¬ 
put  distorts  the  SLL  below  the  main  beam  angle,  60  (i.e.,  in  the  forward  di¬ 
rection),  while  a  flared  output  affects  the  SLL  above  0q  (i.e.,  toward  broad¬ 
side).  By  sacrificing  efficiency  (or  increasing  the  aperture  length),  one  can 
make  the  taper  (or  flare)  more  gradual  to  improve  the  pattern  characteris¬ 
tics.  This  would  be  undesirable  for  most  practical  TWA  structures,  so  we 
consider  only  constant  width  apertures.  Parametric  studies  of  leaky  wave 
guides  with  constant  aperture  width  are  presented  to  develop  a  complete 
design  procedure  for  this  type  of  antenna. 
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Figure  4.  (a)  Perturbed 
waveguide  antenna 
examples:  TE  long  slot 
and  holey  rectangular 
waveguide  showing 
tapered  transition 
sections,  and  (b)  an 
aperture  with  a  constant 
width  section  and  either 
a  tapered  input  or  a 
flared  output  section. 
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Figure  5.  Radiation 
pattern  of  a  leaky 
waveguide  having  an 
aperture  with  tapered 
input  or  flared  output 
sections. 


3.  Perturbation  Solution 


Consider  an  aperture  fed  by  a  parallel  plate  waveguide  radiating  into  a 
half-space  (i.e.,  an  infinite  flange  approximation),  as  shown  in  figure  6, 
where  the  origin  is  centered  in  the  aperture.  We  consider  radiation  into  a 
half-space  rather  than  free  space  in  anticipation  of  constructing  an  antenna 
array  of  leaky  waveguides.  In  the  direction  transverse  to  the  propagation 
direction  y,  the  medium  parameters  over  each  waveguide  cross  section  are 
constant  so  that  a  modal  decomposition  is  possible.  However,  the  wave¬ 
guide  fields  in  the  presence  of  the  aperture  cannot  be  described  by  a  single 
TE  or  transverse  magnetic  (TM)-mode.  We  construct  an  alternate  modal  de¬ 
composition  to  have  the  correct  spatial  variation  and  satisfy  the  boundary 
conditions.  Following  the  analysis  of  Goldstone  and  Oliner  [5],  the  result  is 
a  predominantly  TE-  or  H- type  mode,  as  shown  in  appendix  A.  For  now, 
consider  an  air-filled  waveguide  having  free-space  permittivity,  e  =  eo  and 
permeability,  y  =  po,  propagating  time-harmonic  electromagnetic  (EM) 
fields  with  the  time  dependence,  suppressed.  Since  there  is  no  variation 
of  the  waveguide  in  the  direction  of  propagation,  the  wavenumber  in  this 
direction  is  the  same  for  all  modes  or  ky„  =  ftn  =  ft.  kt o  =  kz o  =  y/k"2  —  ft'2 
is  the  wavenumber  in  the  transverse  (to  y)  direction,  and  Zho  =  uyo/kzo 
is  the  if -type  mode  characteristic  (wave)  impedance.  We  develop  a  com¬ 
posite  equivalent  circuit  to  approximate  the  aperture  of  infinite  length  as 
shown  in  appendix  B.  We  use  the  transverse  resonance  method  described 
in  appendix  C  to  determine  the  transverse  wavenumber  for  the  transverse 
transmission  line  terminated  in  this  equivalent  circuit. 


Figure  6.  A  rectangular 
waveguide-fed  aperture 
radiating  into  a 
half-space. 
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The  transverse  resonance  solution  provides  the  perturbed  waveguide  phase 
constant  with  the  leaky  wave  attenuation  function  (see  sect.  2)  to  com¬ 
pletely  characterize  the  aperture  source  so  that  antenna  gain  can  be  cal¬ 
culated.  We  then  extend  the  result  to  the  finite  length  aperture  by  using  the 
FEM  to  calculate  the  aperture  radiation  efficiency  (see  app  D). 

3.1  Parallel  Plate  Fed  Aperture 

A  parallel  plate  transverse  electromagnetic  (TEM)-mode  polarized  along  x 
is  incident  on  an  infinitely  long  aperture  of  constant  width  radiating  into 
a  half-space  as  shown  in  figure  7.  In  appendix  A,  we  develop  an  alternate 
modal  decomposition  that  can  represent  the  aperture  perturbation  on  the 
waveguide  fields  for  if -mode  or  TE-excitation.  In  appendix  B,  we  show 
that  energy  propagation  in  the  parallel  plate  region  is  equivalent  to  an 
obliquely  incident  TEM  wave  on  the  aperture  plane.  An  if -type  modal  de¬ 
composition  is  used  when  the  TEM  fields  are  the  modal  fields  transverse  to 
the  direction  of  incidence  (see  fig.  B-2).  We  derive  a  variational  expression, 
equation  (B-38)  for  the  radiation  admittance  of  the  parallel  plate  waveguide 
by  assuming  a  constant  aperture  field.  We  then  extend  this  expression  to 
the  case  of  a  narrow  aperture  (i.e.,  when  the  aperture  width  is  not  the  full 
waveguide  height).  The  variational  expression  is  normalized  to  the  char¬ 
acteristic  admittance  of  the  propagating  wave,  which  for  an  incident  Hio- 
mode  is  that  of  the  fundamental  if -type  mode,  so  from  equation  (A-16), 

Ym  =  jt-  =  —  -  (13) 

Zno  w/l  0 


Figure  7.  A  parallel 
plate  waveguide 
radiating  into  a 
half-space. 
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The  variational  result  equation  (B-29)  derived  in  appendix  B  depends  on 
the  E-field  in  the  aperture.  For  a  constant  aperture  field,  one  obtains  the 
approximation  equation  (B-30)  to  the  variational  expression.  This  result 
was  obtained  by  Marcuvitz  using  the  variational  method  [3]  and  includes 
only  the  susceptance  outside  the  aperture,  because  it  depends  only  on  the 
external  E-field.  Then  under  the  same  assumptions,  equation  (B-30)  pro¬ 
vides  the  external  admittance  of  a  narrow  aperture  by  replacing  b  with  d. 
This  is  because  the  power  per  unit  length  radiated  by  the  aperture  is  in¬ 
dependent  of  the  height  of  the  feeding  waveguide,  as  long  as  the  aper¬ 
ture  field  is  unchanged  [1].  For  a  constant  aperture  E-field,  the  normalized 
input  admittance  of  a  parallel  plate  fed  aperture  of  width  d  <  b  is  from 
appendix  B: 


ext  _  _^n_ 

V  m  ~Ym 


(hod)  + 


7rkz()d 


in  which  is  the  Hankel  function  of  the  second  kind,  order  v.  We  use 
the  small  argument  form  of  the  cylindrical  Bessel  functions  [8,  p  935]  in 
equation  (B-38)  to  obtain  a  closed  form  approximation  to  the  external  con¬ 
ductance  and  susceptance  [3,  p  184] 


G,ext 

__  hod 

and 

(15) 

Yho 

Bext 

2  ’ 

Tve 

(16) 

Yho 

7 r 

. gkzod _ 

The  irrational  number,  e  =  2.71828...,  is  the  base  of  the  natural  logarithm 
and  g  =  ec  =  1.781...,  in  which  the  irrational  number  c  =  0.577...,  is  Euler's 
constant.  This  represents  the  "external"  equivalent  circuit  parameters  of  a 
waveguide-fed,  narrow  aperture  radiating  into  a  half-space. 

We  obtain  the  "internal"  circuit  by  considering  the  aperture  as  a  symmet¬ 
rical  capacitive  iris  oriented  along  the  y- axis,  as  shown  in  figure  8a.  The 
internal  fringing  field  of  the  aperture  is  accounted  for  by  the  internal  sus¬ 
ceptance  of  this  symmetrical,  zero-thickness  "window."  The  approach  is  to 


Figure  8.  A  composite 
circuit  model  for 
parallel  plate  fed 
aperture  developed  by 
analog  to  a  capacitive 
iris:  (a)  a  symmetrical, 
capacitive  iris  in 
rectangular  waveguide, 
a  cross-sectional  view, 

(b)  fringing  field  of  an 
iris,  side  view,  and 

(c)  internal  fringing 
field  of  a  parallel  plate 
fed  aperture. 
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consider  the  fringing  fields  on  each  side  of  the  perforated  surface  as  in¬ 
dicated  schematically  in  figure  8b.  The  equivalent  circuit  parameters  for  a 
zero-thickness  window  are  obtained  by  the  equivalent  static  method  for  a 
static  aperture  field  attributable  to  incidence  of  the  two  lowest  //-modes 
[2,  p  219].  The  result  is  the  capacitive  susceptance  given  in  equation  (2a)  of 
The  Waveguide  Handbook,  in  terms  of  the  guide  wavelength  in  the  direction 
of  propagation  [2].  For  our  case  of  a  narrow  aperture,  we  have  d  <C  b,  so 
in  this  expression,  the  second  term  reduces  to  the  function  Q 2  and  the  last 
term  vanishes.  Furthermore,  the  waveguide  height  is  small  compared  to 
wavelength,  or  b  <C  X9  so  that  Q2  ~  0  and  the  second  term  is  negligible. 
Then  for  the  transverse  transmission  line  with  guide  wavelength,  Xgz,  the 
normalized  iris  susceptance  is  the  leading  term  of  equation  (2a)  [2,  p  218] 


-Siris  4  b 

Ym~  V  n 


(17) 


By  symmetry,  we  use  one-half  of  this  result  for  the  iris  in  the  closed  wave¬ 
guide  to  represent  the  internal  susceptance  of  the  aperture  as  indicated  in 
figure  8c. 

The  internal  susceptance  for  the  waveguide  fed  aperture  is  then 


Bmt 

y, m 


kzob  ^ 
7 r 


(18) 


reducing  to  a  short  circuit  when  d  =  0.  We  combine  the  admittance  inside 
and  outside  the  waveguide  to  obtain  a  composite  equivalent  circuit  for  a 
transverse  resonance  analysis.  The  circuit  parameters  are  the  combination 
of  equations  (15),  (16),  and  (17),  or 


G  _  _  kzod 

Yho  Yho  2 


(19) 


B  Bext  Bmti 
Yho=Y^  +  Y^ 
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7re 

.gkzod 


.  kz0b 
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7 r 


(20) 


The  aperture  admittance  for  this  composite  equivalent  circuit  model, 
YA(kzo)  =  G  +  jB,  then  depends  on  the  transverse  wavenumber  and  the 
aperture  geometry.  This  approximation  is  only  appropriate  for  the  ideal¬ 
ized  case  of  zero  wall  thickness  and  an  infinite  waveguide  flange  but  pro¬ 
vides  useful  results  for  narrow  apertures.  A  refined  equivalent  circuit  to 
account  for  the  waveguide  wall  thickness  could  be  developed  but  is  not  re¬ 
quired  at  this  point  to  obtain  engineering  results.  Flowever,  the  equivalent 
circuit  must  be  modified  for  an  array  of  such  leaky  waveguide  antennas, 
including  mutual  interactions  between  nearby  radiating  elements  [4]. 

The  cross  section  of  the  physical  structure  shown  in  figure  9a  is  analyzed 
with  the  transverse  resonance  technique  that  is  described  in  appendix 
C.  We  use  the  network  representation  shown  in  figure  9b  to  terminate 
the  transverse  transmission  line  that  propagates  the  fundamental  H- type 
mode.  A  transverse  field  modal  decomposition  in  terms  of  rectangular 
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Figure  9.  An  infinitely 
long  aperture  in  the 
narrow  wall  of 
rectangular  waveguide 
radiating  into  a 
half-space:  (a)  a 
waveguide-fed  aperture 
and  (b)  equivalent 
circuit  model  for 
transverse  transmission 
line. 


(a)  ,  (b) 
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waveguide  TE-  and  TM-modes  would  require  a  transmission  line  repre¬ 
sentation  for  each  mode  that  is  then  coupled  by  the  aperture.  An  advan¬ 
tage  of  the  if -type  modal  representation  is  that  the  transverse  transmission 
line  consists  of  one  transmission  line  of  characteristic  admittance  Yho,  and 
propagation  constant  jkz o.  This  transmission  line  is  shorted  at  one  end  and 
terminated  in  a  lumped  admittance  at  the  other  end.  The  transverse  reso¬ 
nance  condition  requires  that  at  resonance,  the  input  admittance  "looking" 
in  each  direction  must  sum  to  zero  at  any  point  along  the  line.  We  use  the 
coordinate  system  shown  in  figure  9a,  so  for  any  point  0  <  z'  <  a  we  have 

Ym(z')  +  Yin(z')  =  0,  (21) 

in  which  the  arrows  denote  the  input  admittance  looking  to  the  left  or  right, 
respectively  [7,  p  167]. 

At  the  aperture  plane  (z1  =  a),  the  admittance  looking  to  the  right  is  Ya 
while  that  to  the  left  corresponds  to  a  short  circuit  located  a  distance  "a" 
down  the  line  (i.e.,  at  z'  =  0).  The  resulting  network  resonance  equation  at 
this  terminal  plane  is  then  [2,  p  10] 


G  B 

-j  cot (kz0a)  +  —  +  j—  =  0  ,  (22) 

which  generally  requires  numerical  solution.  An  approximate  analysis  is 
obtained  by  the  perturbation  method  since  the  leaky  waveguide  propaga¬ 
tion  constant  can  be  regarded  as  a  perturbation  on  the  propagation  constant 
of  modes  in  the  closed  waveguide.  This  approach  provides  an  analytical 
approximation  for  the  transverse  wavenumber  as  opposed  to  a  numerical 
evaluation  of  equation  (22).  The  result  is  an  analytical  approximation  to 
the  perturbed  fundamental  mode  propagation  constant,  70  =  j/3,  which  is 
useful  in  the  study  of  leaky  wave  antennas. 

The  analysis  assumes  that  the  aperture  perturbs  the  transverse  wavenum¬ 
ber  from  its  value  in  the  closed  guide.  Then  for  Hio-mode  propagation,  the 
transverse  wavenumber  is  not  much  different  than  k  =  7 r/o — the  unper¬ 
turbed  wavenumber  of  the  lowest  mode.  Then  for  a  small  perturbation,  the 
left  side  of  equation  (22)  is  approximated  by  two  terms  in  a  Taylor  series 

expansion  Y (k)  =  Y (kz 0)  +  Y (kz 0),  at  some  frequency  corresponding  to 
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the  wavenumber,  k.  Then  with  Akzo  =  kz 0  -  k,  we  can  write 


Y{k)  =  Y(k)  + 


dYM 

dkzQ 


Akzo  ~  0 , 

kzO—K 


(23) 


since  for  a  small  perturbation  equation  (22)  will  nearly  vanish  at  resonance. 
Substituting  equation  (22)  evaluated  at  k  into  equation  (23),  we  have 


Akzo  —  —  - 


Y(k) 


-j  cot  (na)  +  Y(k) 


j  Y(k)  -  j  cot (ko) 


dYjkz  q) 
dkzo 


kz  o=k. 


aj  esc  2(«a)  + 


dY(kz0) 

dkZQ 


CSC' 


!(«a) 


(24) 


kzo=K 


The  approximation  is  valid  over  a  wide  range  since  the  derivative  term  is 
inversely  proportional  to  the  fourth  power  of  the  admittance  magnitude 
and  so  is  neglected.  With  the  same  argument,  the  equivalent  expression  in 
terms  of  impedance  quantities  is 


Akzo 


j  Z(k)  +  j  tan(/ca) 
a  sec2(na) 


(25) 


For  a  narrow  aperture  (i.e.,  d/b  <  1),  we  take  the  first  approximation 
to  kzo  as  k  corresponding  to  the  unperturbed  waveguide  propagating  the 
Hio-mode.  Then  with  k  =  ir/a,  the  perturbation  to  the  transverse  wavenum¬ 
ber,  Akzo /  is  calculated.  The  approximate  result  for  such  a  perturbation 
about  a  short  circuit  (i.e.,  the  closed  waveguide)  is  then 


Akzo 
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a  a  \G2  +  B2J 
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The  fundamental  mode  transverse  wavenumber  is  then 

Bf  j  a 


7 r 


kzQ  —  K  +  Akzo  — - h 


+  ^- 


a  a  (G'2  +  B'2)  a  (G'2  +  B'2)  ' 


(26) 

=  K  =  7T /a, 

"  ae 
.gd_ 

•  (27) 

(28) 

The  complex  propagation  constant  of  the  fundamental  mode,  70  = 
\A2-fc*o  =  ar + j(3,  can  then  be  calculated  with  the  root  chosen  so  that  the 
attenuation  is  negative  (ar  <  0)  for  propagation  in  the  positive  y-direction, 
or  (3  =  kyo .  The  propagation  constant  reduces  to  that  of  the  closed  wave¬ 
guide  when  d  —  0  since  the  equivalent  network  becomes  a  short  circuit. 
Note  that  the  result  often  found  in  the  literature  for  leaky  wave  structures 
of  this  type  [1,  p  189]  involves  another  approximation  not  used  here.  The 
propagation  constant  is  often  approximated  by  a  binomial  expansion  with 
higher  order  terms  neglected,  or 
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(29) 


This  approximation  holds  when  the  operating  frequency  is  well  above  the 
cutoff  frequency  of  the  unperturbed  mode  [5]  and  so  does  not  provide  the 
correct  behavior  at  low  frequency.  We  prefer  to  calculate  the  complex  prop¬ 
agation  constant  from  equation  (29)  without  approximation  since  operation 
near  cutoff  is  often  of  practical  interest. 


3.2  Aperture  in  a  Rectangular  Waveguide 

Now  that  70  for  the  parallel  plate  fed  aperture  can  be  calculated  with  suf¬ 
ficient  accuracy  as  a  function  of  /o,  consider  the  finite  length  aperture  in 
a  rectangular  waveguide.  We  assume  that  the  waveguide  width  is  suffi¬ 
cient  to  allow  the  parallel  plate  guide  analysis  to  apply  to  the  rectangular 
waveguide.  This  requires  that  higher  order  modes  associated  with  the  aper¬ 
ture  perturbation  be  rapidly  attenuated  so  that  a  single  propagating  mode 
can  be  used  in  the  transverse  resonance  analysis.  Such  a  leaky  waveguide 
supports  H- type  hybrid  modes  constructed  to  have  the  required  spatial 
variation  (see  app  A).  The  approach  is  to  allow  for  a  longitudinal  variation 
in  d(y)  and  thus  variable  70 (y)  over  the  total  aperture  length,  L.  The  input 
and  output  air-filled  waveguide  propagates  a  single  Hjo-mode  with  prop¬ 
agation  constant  (3  =  /3m  =  \A-o  —  (f)2  for  the  fundamental  {m  =  1 ,  n  =  0) 
mode.  The  Hio-mode  characteristic  impedance  is  Zm  =  but  when  per¬ 
turbed  by  the  aperture,  the  characteristic  impedance  becomes  that  of  the 
,  H- type  fundamental  mode,  Z^o  =  7^7.  This  leads  to  a  nonuniform  trans¬ 
mission  line  representation  along  the  aperture  length 

—fv{y)  =  no(y)ZHoI(y)  and  -  -j-I{y)  =  ho{tj)YmV(y) ,  (30) 

dy  dy 

which  is  solved  numerically  with  the  FEM  with  linear  interpolation  func¬ 
tions  [8,  p  33].  The  //-type  mode  impedance  variation  with  y  is  neglected 
since  for  the  fundamental  mode 


tvz"m  = 


UJ/.IQ 


UJ/IQ 


d 


dy  kz()(y)  k‘i0(y)  dy 


kzo(y)  ~  0  , 


(31) 


because  Re(kzo{y))  >  0  is  a  slowly  varying  function  of  position.  Then  the 
transmission  line  voltage  (and  current)  satisfies  the  differential  equation 


~^y  =  0  5  (32) 

in  which  q(y)  =  -({n/o{y)ZHo{y))~1t  and  C (y)  =  lo(y)/Zno(y)-  The  dis¬ 
tributed  discontinuity  attributable  to  the  (possibly  variable)  aperture  width, 
d(y),  then  causes  a  reflection  in  the  input  waveguide,  as  indicated  in  fig¬ 
ure  10,  where  now  the  aperture  extends  from  y  =  0  to  y  =  L.  For  con¬ 
stant  aperture  width,  the  //-type  fundamental  mode  impedance  is  con¬ 
stant.  However,  there  is  an  impedance  mismatch  at  the  terminal  planes 
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where  the  transmission  line  impedance  is  Z\q  so  that  Zh o(0)  =  Zho(L)  = 
Z10.  Transmission  to  the  matched  output  waveguide  is  also  obtained  in  the 
FEM  solution  with  the  remaining  power  radiated  by  the  aperture.  For  a 
single  incident  mode  with  unit  input  voltage,  the  total  voltage  for  y  <  0  is 
then 

V(y)  =  +  Tio ,  (33) 

in  which  T iq  is  the  voltage  reflection  coefficient  at  the  input  terminal  plane 
of  a  finite  length  section  of  leaky  waveguide.  The  boundary  conditions  at 
the  terminal  planes  are  determined  from  the  functional  for  the  differential 
operator  in  equation  (32)  [8,  p  34].  Minimizing  the  functional  establishes 
the  matrix  equation  to  solve  where  the  boundary  terms  can  be  cast  in  the 
form  of  boundary  conditions  of  the  third  kind  (or  mixed  boundary  con¬ 
ditions).  However,  the  derivative  terms  at  the  boundary  y  —  0  are  taken 
in  the  negative  sense  so  that  in  general,  boundary  condition  is  written  as 
±q(dV/dy)  +  yF  =  p.  The  negative  sign  is  used  at  y  =  0  so 


-C 


dV 

dy 


~J  7o 


y= 0  7oZtfo(2/) 


(e~rm  _  Twenoy)j 


V— o 


Zho{v) 


(e~noy  _  rhoe™’3')  ,  (34a) 


y= o 


and  we  can  identify 

Xo  =  and  p0  =  .  (34b) 

z  io  Xio 

The  subscript  denotes  the  input  terminal  plane  at  y  =  0  where  Zho(0)  = 
Ziq.  Similarly,  for  y  >  L,  the  voltage  is  V(y)  =  Twe~^°y  in  which  T10  is  the 
slotted  waveguide  transmission  coefficient.  Then  at  y  =  L,  we  have 
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dV 
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-no 


y=L  -loZmiy) 

Thus,  we  can  identify 


(Tioc“r,oy) 


y=L 


Zho(v) 


Tioe' 


-nov 


y=L 


(35a) 


XL  =  ~~i~  and  pL=0  (35b) 

^10 

at  the  end  of  the  finite  length  aperture  where  Zhq(L)  =  Z\q.  The  trans¬ 
mission  line  voltage,  V (y)  is  obtained  directly  from  the  FEM  solution  once 
7o(y)  is  input  from  the  transverse  resonance  result  for  a  given  aperture 
width  profile,  d{y).  The  distributed  impedance  discontinuity  in  the  leaky 
waveguide  section  reflects  some  power  to  the  input  waveguide  and  trans¬ 
mits  some  power  to  the  output  waveguide  assumed  to  be  terminated  in  a 
matched  load  (or  infinite  in  extent). 

The  transmission  line  model  characterizes  propagation  on  a  nonuniform, 
lossy  transmission  line  by  the  variation  in  70(2/)  over  the  length  of  the 
waveguide  section  (for  variable  d{\j)).  The  power  lost  because  of  dissipa¬ 
tion  in  the  leaky  waveguide  is  considered  completely  radiated.  At  y  =  0, 
the  voltage  is  given  by  equation  (33)  so  that  710  =  V(0)  -  1.  At  y  =  L, 
V(L)  =  Tioe“-?7o^)L  with  70 (L)  =  7o(0)  =  /?io  at  each  end  of  the  leaky 
waveguide  section.  The  reflection  and  transmission  coefficients  are  deter¬ 
mined  from  the  FEM  voltage  solution.  We  define  the  antenna  efficiency  as 
the  combination  of  the  input  efficiency  (reflection  loss),  the  radiation  effi¬ 
ciency  (load  dissipation)  and  other  losses  combined  into  t]c  (typically  dom¬ 
inated  by  conductor  losses) 

V  =  VinVrVc  =  (l  -  |r10|2)  (1  -  |T10|2)  (1  -  |  A10|2)  (36) 

with  the  absorption  coefficient,  A 10,  presently  neglected  or  r)c  =  1-  We  in¬ 
vestigated  the  convergence  of  the  FEM  solution  by  increasing  the  number 
of  cells  used  per  guide  wavelength.  The  results  showed  that  a  FEM  mesh 
with  more  than  75  cells  per  guide  wavelength  is  required  for  the  solution 
to  converge.  For  the  results  presented,  we  used  an  FEM  mesh  with  100  cells 
per  guide  wavelength  to  calculate  y  (see  app  D).  A  tapered  transition  region 
to  the  constant  width  section  does  reduce  F 10,  but  the  effect  is  insignificant 
for  a  narrow  aperture.  For  the  finite  length  apertures  considered,  the  effect 
of  tapered  input  and  output  transition  sections  on  Tjo  is  negligible  but  can 
distort  the  radiation  pattern  as  previously  mentioned. 
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4.  Radiation  Pattern 


The  radiated  fields  produced  by  arbitrary  current  distributions  can  be  cal¬ 
culated  with  the  appropriate  Green's  function,  which  is  the  radiated  field  of 
an  infinitesimal  source  [8,  p  909].  For  a  rectangular  aperture  symmetrically 
located  in  the  coordinate  system  of  figure  6  (also  see  fig.  1),  the  far  fields  in 
free  space  are  computed  with 


2ie-jfcrx 


f  j  f2d  (cos 9 sin 0P  -  cos 4>t)Ex(x\y')e^ksines'in^3(y'^y' ejksin0cos^x' dx'dyf 


,A.j  (  >  and 


o  -jkr  y 
4tt  r  e  X 


/f  f2d  (cos 9 sin <pt  +  cos<frp)Ex(x\y’)ertksinesin*-P{-y,')'>v' eiksmdcos*x' dx'dy' 


where  t  and  p  are  unit  vectors  in  the  9  and  4>  directions,  respectively.  Here 
Ex(x' .  y')  is  the  aperture  E-field  and  /3(y')  is  the  perturbed  waveguide  phase 
constant  at  each  position  y' .  In  the  far  zone,  the  EM  fields  are  of  the  form  ' 


■VIoE(9,(t>)xT  and  H(r) 


M),  (38) 


in  which  =  120-7r-ohm  is  the  free-space  impedance  and  r/o  =  l/<fo  is  the 
free-space  admittance.  The  average  power  radiated  per  unit  area  is  then 

S  =  ExH<=(FxtrXF=W2,  (39) 


so  |F(0^)|2  is  the  average  power  radiated  per  steradian.  The  directive  gain 
is  the  radiated  power  per  steradian  normalized  to  the  total  radiated  power, 

f  r  r2ir  rTT  2 

Pt otai=  //  IF  (0)1 2dn=  /  /  \F(6,<f>)\  sin 9d9d4>  (40) 

JJiit  J  0  Jo 

averaged  over  all  solid  angles  (i.e.,  47t  steradian).  Therefore,  the  directivity 
(or  directive  gain)  is  given  by 


G{9,<j>)  = 


4tt|  F{9,(p): 


in  dB  relative  to  an  isotropic  radiator  (dBi)  and  is  independent  of  amplitude 
scale  factors.  The  realized  gain  (or  simply  gain)  includes  the  antenna  effi¬ 
ciency  as  calculated  with  the  FEM  (see  sect.  3.2).  The  realized  gain,  Gr(9,  </>), 
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includes  losses  such  as  the  input  impedance  mismatch  (?/ jn),  the  power  dis¬ 
sipated  in  the  load  (»/,.),  and  other  possible  loss  mechanisms  such  as  con¬ 
ductor  or  dielectric  loss  (i]c).  The  calculated  results  are  presented  as  realized 
gain,  based  on  the  total  antenna  efficiency,  or  Gr(9, 4>)  =  r/G(9,  </>). 

The  finite  width  rectangular  aperture  can  be  represented  by  an  equivalent 
rectangular  sheet  of  current  extending  from  x  —  -d/2  to  x  =  d/2  and 
y  =  —L/2  to  y  =  L/2  (see  fig.  la).  Since  there  are  no  transverse  currents 
excited  for  the  assumed  aperture  .E-field,  the  amplitude  depends  only  on 
axial  position,  Ex(x',y')  =  Er(y').  Then  even  for  variable  aperture  width, 
d(y'),  the  integral  in  equation  (37)  along  the  .r-direction  simplifies  to 


.Odl  4G1 

12a  EA^'.y'y^'^'dx'  =  e,W) 


sin  6 cos 4IT' £  f  /s£^(^x(y  )) 

*  J  ux{y') 


(42) 


in  which  ux(y')  —  it d(y')  sin  9  cos  (t)/\n.  However,  we  only  consider  narrow 
apertures  in  which  d(y')  <  A0/2  so  ux(y')  =  ux  <  1  and  the  aperture 
field  in  equation  (37)  reduces  to  a  magnetic  current  line  source  as  shown  in 
figure  11  where 

M  =  -z  x  E  =  -Ex(y') y  .  (43) 


Then  in  equation  (42),  the  radiation  pattern  is  not  modified  by  the  finite 
aperture  width  as  long  as  it  is  small  compared  to  wavelength  (i.e.,  d  <  Ao/2) 
[1,  p  47], 

For  constant  aperture  width,  ux{y')  =  ux  =  nd  sin  9  cos  <£/%,  and  Ex(x',  y')  = 
Ex{y')  The  line  source  amplitude  is  the  aperture  illumination  function  mag¬ 
nitude  from  equation  (12),  or 

d 

Er(y')  =  |f(j/)|  [l  ejks>ne<:os&r' dx'  =  |f(?/)|  ^  =  |f(?/)|  (44) 

J~  2  Ux 
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This  aperture  field  has  the  phase  variation  of  the  perturbed  waveguide 
phase  constant  which  is  now  independent  of  position  in  the  aperture, 
/ 3(y ')  =  (3.  Here,  we  consider  constant  width  apertures  so  that  results  for 
variable  aperture  width  are  shown  only  in  figure  5.  For  the  narrow  widths 
of  interest,  the  far  fields  in  the  upper  half-space  are  from  equation  (37) 

E(r)  «  {  iLe-jkr  r  (cos  0  sin  <f>ft  -  cos  <j>t)Ex (y  V(fc sin 6 sin dy'  1  ,  and  (45a) 
{2irr  Jyo  J 

H(r)  »  jj  (cos  0  sin <f>t  +  cos  4>p)Ex (y')eJ(ksm 9 ^~P)y' dy'  j  .  (45b) 

These  results  are  sufficient  for  radiation  pattern  calculations  since  the  di¬ 
rectivity  from  equation  (41)  is  independent  of  amplitude  scale  factors. 

We  conduct  a  parametrical  study  to  evaluate  the  influence  of  slot  length 
and  width  on  the  far  field  radiation  pattern  at  fixed  frequency.  For  the 
waveguide  and  slot  dimensions  previously  considered  ( ka  =  4.25,  b  =  a/2 
and  d  =  b/ 30),  the  H- plane  (i.e.,  4>  =  7r/2)  realized  gain  pattern  is  shown 
in  figure  12  with  aperture  length  as  a  parameter.  The  total  efficiency  is  also 
shown  corresponding  to  rj  —  —1.4,  —0.8,  —0.5,  and  —0.3  dB.  The  pattern 
is  as  expected  for  an  approximately  uniform  source  distribution  with  SLL 
~  13  dB.  The  beam  appears  at  an  angle  that  depends  on  the  phase  varia¬ 
tion  over  the  aperture,  and  this  is  completely  determined  from  the  variable 
aperture  width.  The  azimuthal  beam  angle,  4>o  =  7t/2,  remains  fixed  and 
is  referred  to  as  the  "forward"  direction.  For  constant  aperture  width,  the 
elevation  beam  angle,  6o,  is  determined  by  the  perturbed  waveguide  wave¬ 
length.  As  the  aperture  width  approaches  zero,  the  elevation  beam  angle 
approaches 

60  =  sin”1  (%/Aa0)  ~  42°  (46) 

in  which  Xgo  is  the  guide  wavelength  for  the  unperturbed  (i.ev  closed) 
waveguide.  This  is  shown  in  figure  13  for  fixed  slot  length,  L  =  20Ao,  and 


Figure  12.  Realized  gain 
for  a  constant  width 
aperture  in  a  waveguide 
narrow  wall  with  length 
as  a  parameter. 
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Figure  13.  Realized  gain 
for  a  constant  width 
aperture  in  a  waveguide 
narrow  wall  with  width 
as  a  parameter. 
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decreasing  slot  width  where  the  realized  gain  is  reduced  according  to  the 
total  efficiency. 

Now  consider  a  specific  example  of  a  constant  width  aperture  in  the  narrow 
wall  of  WR-284  waveguide  where  the  interior  dimensions  are  a  =  2.84  in. 
and  6  =  1.34  in.  (which  is  not  half-height  waveguide).  The  aperture  width 
is  constant  at  d  —  b/ 25  with  length  L  =  30o  ~  20Ao  at  a  frequency  where 
ka  =  4.25.  The  calculated  gain  at  this  frequency  is  shown  as  a  function  of 
the  elevation  angle  and  versus  the  azimuthal  angle  at  fixed  elevation.  The 
H -plane  gain  with  <f> o  =  7r/2  is  shown  in  figure  14a  where  the  main  beam 
angle,  9q  =  24.7°. 

The  realized  gain  versus  azimuth  at  this  fixed  elevation  angle  is  shown  in 
figure  14b.  The  beam  can  be  frequency  scanned  from  broadside  and  re¬ 
mains  free  of  grating  lobes.  The  main  beam  angle  from  broadside  (or  over¬ 
head)  when  4>o  =  n/2  is  shown  in  figure  15  as  a  function  of  normalized  fre¬ 
quency.  The  realized  gain  for  WR-284  slotted  waveguide  (with  d  =  6/30)  is 
shown  in  figure  16  at  three  frequencies  corresponding  to  ka  =  4.0, 4.25,  and 
4.5,  in  which  t]  =  —0.9.  —1.3,  and  —1.7  dB,  respectively.  The  result  indicates 
that  the  pattern  is  stable  over  this  frequency  BW  with  less  than  a  2-dB  gain 
variation.  The  main  beam  can  be  positioned  at  a  desired  elevation  angle  by 
adjusting  the  guide  phase  constant  and  then  frequency  scanning  over  about 
20°  (for  12  percent  BW).  The  corresponding  peak  gain,  Gr(9o,  n/2)  and  total 
efficiency  are  shown  in  figure  17a  and  17b,  respectively.  To  obtain  the  radi¬ 
ation  pattern  for  a  realistic  aperture,  the  rectangular  source  distribution  in 
two  dimensions  is  required  as  in  equation  (37).  However,  for  the  narrow 
slots  considered  (without  transverse  currents),  equation  (45)  is  used  since 
the  radiation  by  a  rectangular  aperture  source  is  a  negligible  correction  of 
these  line  source  results. 
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Figure  14.  (a)  iJ-plane 
realized  gain  for  a 
constant  width  aperture 
in  the  narrow  wall  of 
WR-284  and  (b)  realized 
gain  at  the  main  beam 
angle  for  a  constant 
width  aperture  in 
narrow  wall  of  WR-284. 


Elevation  angle  0° 


(b) 
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Figure  15.  Main  beam 
angle  scan  range  for  a 
constant  width  aperture 
in  narrow  wall  of 
WR-284. 


Figure  16.  Calculated 
H -plane  gain  for  an 
aperture  in  narrow  wall 
of  WR-284  waveguide. 
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Figure  17.  Results  for  a 
constant  width  aperture 
in  narrow  wall  of 
WR-284  waveguide:  (a) 
peak  realized  gain  and 
(b)  total  efficiency. 


(a) 


5.  Conclusion 


Following  Goldstone  and  Oliner  [5],  we  have  presented  an  equivalent  cir¬ 
cuit  model  for  a  narrow  slot  in  a  rectangular  waveguide.  The  transverse 
resonance  technique  is  applied  to  the  resulting  transmission  line  descrip¬ 
tion  to  provide  the  resonance  condition.  The  transcendental  equation  is 
then  approximated  by  a  perturbation  solution  resulting  in  the  perturbed 
waveguide  propagation  constant.  The  waveguide  attenuation  constant  de¬ 
termines  the  aperture  distribution  function  that  has  the  waveguide  longi¬ 
tudinal  phase  variation.  The  radiated  EM  fields  are  then  obtained  by  nu¬ 
merical  integration  over  the  aperture  source,  where  the  main  beam  angle 
depends  on  the  perturbed  waveguide  phase  constant.  It  was  shown  that  a 
tapered  (or  flared)  slot  width  corrupted  the  radiation  pattern,  so  we  con¬ 
sider  only  constant  width  slots  in  detail.  Parametric  results  were  shown  for 
a  constant  width  slot  in  WR-284  waveguide,  which  corresponds  to  a  ta¬ 
pered  aperture  illumination  function.  The  TE-slotted  waveguide  TWA  has 
good  radiation  characteristics  (i.e.,  ?/  >  70  percent)  over  a  frequency  range 
corresponding  to  ko  =  4.0  —  4.5  or  rather,  ka  =  4.25  ±  6  percent.  Over  this 
12-percent  frequency  band  (i.e.,  2.65  to  2.98  GHz),  the  main  beam  scans  10° 
to  30°  from  broadside  (or  overhead).  This  scan  range  could  be  augmented 
with  electronic  control  of  the  waveguide  phase  constant  in  order  to  phase 
scan  the  beam  from  broadside.  In  the  sequel  to  (i.e.,  Part  II)  this  report,  di¬ 
electric  loaded  waveguides  are  considered  in  which  voltage-controlled  di¬ 
electric  materials  integrated  into  the  antenna  structure  may  allow  efficient 
phase  scanning  of  the  slotted  waveguide  TWA. 
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Appendix  A.  Modal  Decomposition  in  a  Parallel  Plate  Region 


Consider  a  parallel  plate  region  of  infinite  extent  with  the  coordinate  sys¬ 
tem  shown  in  figure  A-l,  where  we  presently  locate  the  origin  on  the  bot¬ 
tom  plate.  We  desire  a  modal  decomposition  for  electromagnetic  (EM)  fields 
having  an  exponential  phase  variation,  e~^nV  in  the  propagation  direction 
for  the  nth-mode.  A  modal  decomposition  in  terms  of  transverse  to  y  elec¬ 
tric  (TE-)  and  magnetic  (TM-)  fields  is  possible,  but  ultimately,  we  desire 
to  include  an  aperture  in  the  waveguide  which  perturbs  the  modal  fields. 
In  a  typical  TE-TM  mode  decomposition,  separate  transmission  lines  in  the 
direction  transverse  to  y  are  required  to  represent  each  mode,  which  are 
coupled  by  the  presence  of  the  aperture.  We  indicate  this  schematically  in 
figure  A-2,  where  an  aperture  fed  by  a  parallel  plate  waveguide  creates 
fringing  fields  that  have  Ex-  and  f?2 -components.  The  mode  would  also 
have  an  Cy-component  requiring  a  linear  combination  of  TE-  (or  fT-modes) 
and  TM-modes  to  represent  the  EM  field.  Single  mode  decomposition  is 
more  attractive  in  order  to  obtain  a  single  transverse  transmission  line  rep¬ 
resentation  for  an  aperture  in  rectangular  waveguide.  To  this  end,  we  use 
an  alternate  modal  decomposition  that  will  be  helpful  in  the  analysis  of 
leaky  wave  guiding  structures  such  as  that  shown  in  figure  A-3. 

We  desire  to  construct  a  transmission  line  description  of  waves  propagating 
in  a  direction  transverse  to  the  propagation  direction,  y  in  the  parallel  plate 
guide  (i.e.,  the  ^-direction  operator  Vt2  =  From  the  Maxwell 

field  equations,  the  EM  fields  must  satisfy 


dE  u 
dz 


=  juy 


Y  Vt2Vt2 
+  k2 


■  (Ht2  x  z) 


(A-la) 


dH u 
dz 


=  joJE 


j  .  V  tzS? tz 


•  (z  x  E tz) 
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Figure  A-2.  Parallel 
plate  fed  aperture. 


Figure  A-3.  An  aperture 
in  narrow  wall  of 
rectangular  waveguide. 


in  which  I  =  xx  +  yy  +  zz  is  the  unit  dyadic  and 

juj/iH -  =  V/2  •  (z  x  Ef-)  (A-2a) 

jueEz  =  Vtz  ■  (H^  x  z)  .  (A-2b) 

In  the  parallel  plate  region,  0  <  x  <  b,  the  medium  has  permeability,  /i  = 
Po/'r  and  permittivity,  e  =  (oer,  which  are  presently  taken  to  be  constants. 

In  the  spirit  of  separation  of  variables,  expand  the  transverse  (to  z)  fields 
in  vector  modal  functions  that  are  constructed  to  have  the  correct  x  and  y 
spatial  variation. 
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E  tz{x,y,z)  =  Vn(z)en(x,y)  (A-3a) 

H tz(x,  y,  z )  =  In(z)hn(x,  y)  .  (A-3b) 

The  Vn{z)  and  In(z)  are  the  nth-mode  amplitude  functions  that  satisfy  trans¬ 
mission  line  equations  in  the  longitudinal  direction, 

-  —j—  =  jkznZnIn(z)  (A-4a) 

_dl^z)_  =  jkznYnVn{z)  (A-4b) 

dz 

in  which  Zn  =  1/Yn  is  the  characteristic  (wave)  impedance  and  jkzn  the 
propagation  constant  along  z  corresponding  to  the  nth-mode.  When  we 
use  equation  (A-3)  in  equation  (A-l)  and  take  account  of  equation  (A-4), 
the  transverse  (to  z)  mode  functions  are 
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UOfl 

kzn^n 


Y  Vt2Vt2 
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(hn  x  z) 
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(A-5b) 


The  transverse  EM  fields  are  then  obtained  from  equation  (A-3)  in  terms  of 
the  transmission  line  voltage  and  current  functions.  With  equation  (A-3)  in 
equation  (A-2),  the  "longitudinal"  fields  are 


hn  — 


=  Vn{z)Vu  •  (z  x  en)  (A-6a) 

jueEz  —  In[z)^ tz  ■  (hn  x  z)  .  (A-6b) 

With  en  =  h„  x  z,  the  classical  TE-TM  modal  decomposition  could  be  ob¬ 
tained.  Instead,  we  define  the  en-mode  function  by 


-i/?  ii  ( nnx\  „  ,  „ 

en  =  e  3Pny  cos  (  — —  )  x  =  if?nx  . 


(A  -7) 
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Substituting  equation  (A-7)  into  the  right  side  of  equation  (A-5b),  we  obtain 
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With  k 2  =  to2fie  and  \bn  =  e  3/3nV  cos  (^j2),  this  becomes 
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■  (A-8) 
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The  normal  component  vanishes  at  x  =  0  and  x  =  b,  as  required  by  a 
perfectly  conducting  boundary.  To  show  that  equations  (A-9)  and  (A-7)  are 
mutually  consistent,  substitute  equation  (A-9)  into  the  right  side  of  equa¬ 
tion  (A-5a)  to  obtain 
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When  we  take  derivatives  and  cancel  terms,  equation  (A-10)  becomes 


UJfl 
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Then  the  result  in  equation  (A-ll)  must  equal  the  left  side  of 
equation  (A-5a).  This  gives  the  condition 

A *n  =  k*-fi-(jf)2  (A- 12) 

on  the  wavenumber  for  the  transverse  (to  y)  transmission  line.  The  nor¬ 
malization  factor  is  the  characteristic  impedance,  chosen  so  that  the  y- 
component  of  the  h„  mode  function  is  simply  -i pn.  From  equation  (A-9), 
this  requires  that 


k2  -  d2 


=  1  or  Zn  — 


hzrtUJ 


(A-13) 


With  this  normalization,  the  electric  mode  functions  are  from  equation  (A-7) 

(A-14a) 

and  from  equation  (A-9),  the  magnetic  vector  mode  functions  are 


h 


n 


(A-14b) 


The  transverse  fields  are  found  from  equation  (A-3)  and  are  simply  the 
transverse  mode  functions  with  transmission  line  amplitudes.  The  corre¬ 
sponding  longitudinal  components  are  obtained  when  equation  (A-14)  is 
used  in  equation  (A-6)  and  when  we  take  account  of  equation  (A-4).  The 
result  is 


Hz„  =  -  vn(z)  cos  (^pj  e-^y  (A-15a) 
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This  mode  is  a  linear  combination  of  a  TE-mode  and  a  TM-mode,  often 
referred  to  as  an  AA-type  hybrid  mode.  The  AA-type  mode  characteristic 


impedance  is  denoted  Z#n,  to  have  a  clear  distinction  from  the  TE  (or 
il-mode)  case.  We  find  the  fundamental  mode  impedance  from  equation 
(A-13)  using  equation  (A-12)  with  n  =  0, 


Zho  = 


k'lcsUJU  LOU 


(A-16) 


When  n  =  0  in  equations  (A-14)  and  (A-15),  the  H- type  mode  functions 
degenerate  into 


eo  =  xe  ■?/3oy  ,  and 

(A-17a) 

h0  =  ye-jpoy  ,  with 

(A-17b) 

Hon  =  ~  V0(z)e-W°y 

(A-17c) 

Thus,  the  fundamental  H- type  mode  is  equivalent  to  an  Hio-mode  propa¬ 
gating  in  the  y-direction.  The  H- type  mode  decomposition  is  useful  in  the 
analysis  of  leaky  wave  antennas  such  as  an  aperture  in  the  narrow  wall  of 
rectangular  waveguide  (see  fig.  A-3).  The  H- type  mode  decomposition  in 
the  transverse  direction  has  a  fundamental  mode  that  can  represent  an  inci¬ 
dent  Hio-mode.  The  aperture  perturbation  on  this  incident  mode  can  then 
be  described  as  introducing  higher  order  H- type  modes. 
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Appendix  B.  Equivalent  Circuit  Parameters  for  a  Leaky  Hi0 
Rectangular  Waveguide 


We  desire  equivalent  circuit  parameters  for  a  parallel  plate  waveguide  fed 
aperture  radiating  into  a  half-space  as  shown  in  figure  B-l.  To  follow  con¬ 
vention,  the  origin  is  presently  located  on  the  bottom  plate.  The  case  of 
interest  is  an  aperture  in  the  narrow  wall  of  a  rectangular  waveguide  when 
excited  by  the  Hio-  (or  TE10-)  mode.  The  unperturbed  waveguide  prop¬ 
agates  the  Hio-mode,  but  when  the  waveguide  is  perturbed  by  an  aper¬ 
ture,  an  JT-type  hybrid  mode  is  appropriate  as  developed  in  appendix  A. 
The  H- type  modal  decomposition  in  the  transverse  direction  reduces  to  the 
Hio  case  when  n  =  0  as  shown  in  equation  (A-17).  Then,  TE-excitation  can 
be  viewed  as  a  fundamental  ff-type  mode  incident  on  the  aperture  plane, 
and  the  aperture  perturbation  produces  higher  order  modes.  Energy  prop¬ 
agation  along  y  in  the  parallel  plate  region  can  be  viewed  as  a  TEM  wave 
incident  at  some  angle,  6  to  the  aperture  ( xy )  plane  as  shown  in  figure  B-2. 

As  in  appendix  A,  take  the  2-axis  as  a  propagation  direction  and  consider 
transverse  field  components  incident  on  the  ary-plane.  The  electric  (E-)  field 
is  polarized  along  the  ar-axis  and,  is  incident  at  an  oblique  angle  to  the  aper¬ 
ture  plane  (see  fig.  B-2) 


Einc  =  ±e~yky0y~jkz0Z 


(B-l) 


Figure  B-l.  TE-excited 
aperture  radiating  into  a 
half-space  (side  view). 
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Figure  B-2.  TEM-wave 
incident  at  an  oblique 
angle. 


From  Faraday's  Law,  the  incident  magnetic  ( H- )  field  is 


-  y  ^2°  e-jkyoy-jk:oz  _  ^  e-jkyoy-jk:(,z 

to  Ho  to  no 


in  which  only  the  transverse  (to  z)  component  is  required  in  this  analysis. 
For  an  air-filled  waveguide,  the  components  of  the  propagation  constant 
are  jky o  =  j(3 o  =  jk  cos  6  and  jkz0  =  jk  sin  9.  The  wavenumber  magnitude 
is  then 


with  longitudinal  wavenumber,  kz o,  of  the  transverse  (to  y)  transmission 
line. 


These  incident  electromagnetic  (EM)  fields  correspond  to  the  fundamental 
(n  =  0)  mode  of  an  H- type  modal  decomposition  as  developed  in  appendix 
A.  Using  Emc  from  equation  (B-l)  in  equations  (A-lb)  and  (A-2b),  the  inci¬ 
dent  magnetic  ( H- )  field  is  that  of  the  fundamental  FT-type  mode  equivalent 
to  equations  (A-17b)  and  (A-17c),  but  propagating  at  some  oblique  angle  to 
the  z-axis.  This  is  because  the  mode  amplitudes  satisfy  the  transmission  line 
equation  (A-4)  so  the  fields  that  are  transverse  to  the  propagation  direction 
are  related  by  the  //-type  mode  characteristic  (wave)  impedance  equation 
(A-13).  Then  the  transverse  component  in  equation  (B-2)  can  be  written  in 
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terms  of  the  fundamental  mode  wave  impedance,  Zho  =  wno/kzo.  Higher 
order  E- type  modes  for  the  variation  along  x  also  satisfy  the  boundary  con¬ 
ditions  on  the  perfectly  conducting  surfaces  when  kxn  =  (nir/b).  Since  there 
is  no  variation  of  the  waveguide  along  y,  the  wavenumber  in  this  direction 
is  the  same  for  all  n  or  pn  =  /3o  =  /?.  For  the  nth-mode,  the  wavenumber 
magnitude  of  the  transmission  line  along  the  2-axis  is  thus 

fc«,  =  y (B-4) 


The  transverse  modal  decomposition  in  appendix  A  showed  in  equation 
(A-17)  that  the  fundamental  E- type  mode  degenerates  to  an  H10-mode 
propagating  along  y.  Then  the  incident  EM  fields  that  are  transverse  to  2 
can  be  described  by  an  Hio-mode  propagating  along  kinc  (see  fig.  B-2).  The 
incident  EM  fields  transverse  to  2  can  be  written 


Einc  = 


and 


(B-5a) 


P-jpy-jkz0z 

H,nc  =  y - -= -  (B-5b) 

6m 

corresponding  to  a  TEM  wave  incident  at  an  oblique  angle.  These  inci¬ 
dent  fields  are  partially  reflected  along  the  aperture  with  higher  order 
modes  generated  according  to  the  aperture  reflection  coefficient,  Tn,  for 
each  E— type  mode.  The  total  EM  fields  within  the  waveguide  are  then  a 
superposition  of  the  incident  field  and  an  infinite  sum  of  reflected  modes. 
The  total  internal  (i.e.,  2  <  0)  E-field  transverse  to  2  is 


Kz  =  X 


(e-jk*°z  +  r0e^°*)  +  rn  COS  (  ^  J  ejk*nZ 

n—  1  '  ' 


o~jPv 


(B-6a) 
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tz 


The  total  transverse  E-field  is  found  with  equation  (A-5b)  for  the  trans¬ 
verse  mode  functions  by  the  inclusion  of  the  sum  of  reflected  modes.  The 
result  is  as  in  equation  (A-14b)  for  the  h„-mode  function  or 


^  (e-jfczCV-roe^o2) 

^  ZH0 

_  rne]^"z 

n=l  Hn  LK0 


/32 


(^)  Sin(^f)x  +  cos(^f)y 


>  e 


-iPy 


(B-6b) 


where  the  wave  impedance  is  that  for  the  E-type  waveguide  mode 

17  _  kznuiyo  _  kznojjj,o 

',=^=wW'  <B'7) 

The  longitudinal  (2-)  field  components  can  be  found  with  equation  (A-6) 
but  are  not  required  here. 

For  2  >  0,  the  total  EM  fields  can  be  written  as  a  superposition  of  complex 
exponential  functions  corresponding  to  a  Fourier  integral  with  the  trans¬ 
form  variable  being  the  transverse  wavenumber.  With  longitudinal  and 
transverse  components  of  the  /.: -vector,  k  =  k2  +  kt,  the  radiated  fields  are 


35 


and 


(B-8a) 


Eext  =  [  [  E  fk;) 

J—O C  J  -OC  (277 ) ^ 

Hext  =  1°°  f00  ^  X  e-jkxx-jkay-jy/^-k^-klz^R^  (B-8b) 
J—oc  J—o o  (27r)2 

The  wavenumber  along  the  aperture  is  that  of  the  waveguide  mode,  ky  —  fi 
for  all  n  or 

E(kt)  =  2irE(kx)6{ky-0)  .  (B-9) 

Then  the  external  EM  fields  are  found  from  equation  (B-8)  via  the  ampli¬ 
tude  coefficients  in  equation  (B-9) 
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where  the  longitudinal  wavenumber  magnitude  is 


Ax.  =  A’2  -  A'2  —  ft2  or  A2  +  A2  =  k2  —  (32  . 


(B-lOa) 

(B-lOb) 

(B-ll) 


Since  the  component  of  the  wavenumber  along  the  aperture  length,  /?,  is 
fixed,  the  amplitude  coefficient  can  be  written 


E(AX)  —  £x ( kx ) x  +  £z(kx)z  . 


(B-12) 


Further,  there  is  no  component  of  jS-field  along  the  A -vector  in  the  radi¬ 
ation  direction,  k,  so  we  can  eliminate  one  component  of  the  amplitude 
coefficient  by  recognizing  that  k  •  E  =  0  or 


k  £x  +  kXz  =  0  so  £z  =  -[  ^  )  £x  . 


k~ 


(B-13) 


The  cross  product  in  equation  (B-8b)  is  then 


k  x  E  =  [£~/?x  +  (k:£x  —  kx£z)  y  —  f3£xz] 
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■—) 


(B-14) 


The  external  EM  fields  can  be  written  in  terms  of  this  single  component  of 
the  Fourier  transform  variable,  with  the  result 
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Then  the  external  fields  that  are  transverse  to  z  are 
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The  boundary  condition  on  the  internal  and  external  fields  at  the  aperture 
plane  requires  continuity  of  the  tangential  (i.e.,  transverse  to  z)  field  compo¬ 
nents.  Then  equating  equations  (B-6)  and  (B-16)  and  omitting  the  common 
factor  e~J/3y,  we  must  have  at  z  =  0 


OO  /  \  ”|  aQQ 

1  +  ro  +  E  rn  cos  (=)  =  ^j  £xe-jkxXdkx  ,  (B-17a) 
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using  the  relation,  k2 0  =  k2  —  (32.  For  the  nth-mode  we  have 


k zn  —  \  k20 


nir 
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(B-17b) 

(B-17c) 


)  -V(t)  -‘2si“2*  <B'18) 


Then  using  equation  (B-18)  along  with  equation  (B-7),  we  find  that  equation 
(B-17b)  becomes 
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sm 
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kz 


(B-19) 


The  aperture  E-field  at  z  =  0,  ExA(x, 0),  is  the  one-dimensional  Fourier 
transform  that  appears  on  the  right  side  of  equation  (B-17a).  We  have  the 
Fourier  transform  pairs 


1  P°° 

ExA(x)  =  —  /  £xe~jkxXdkx  ,  and  (B-20a) 

Z7T  J  —  oo 

£x  =  E  J  ExA{x')e?k*x' dx’  (B20b) 

Ap 

in  which  the  second  integral  extends  only  over  the  aperture.  Using  equa¬ 
tion  (B-7),  we  find  that  the  condition  equation  (B-17c)  becomes 
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in  which  Hi  is  the  Hankel  function  of  the  second  kind,  order  /a  An  expres¬ 
sion  for  r„  can  be  obtained  from  equation  (B-17b),  but  this  is  not  required 
when  one  is  using  orthogonal  mode  functions. 

We  can  take  advantage  of  orthogonal  functions  to  further  simplify  the 
results.  For  now,  assume  that  the  aperture  is  the  full  waveguide  open¬ 
ing  or  d  =  b  in  figure  B-2.  Multiply  both  sides  of  equation  (B-17a)  by 
cos(rmx/b),n  >  1  and  integrate  over  the  aperture  width,  b.  The  result  is 
an  expression  for  T„ 

^r„  =  ExA(x')  cos  dx'.  n  >  1  .  (B-22a) 


Taking  the  complex  conjugate  of  equation  (B-22a),  we  have  also 

=  j'  E*xA{x')  cos  dx n  >  1  .  (B-22b) 

By  directly  integrating  both  sides  of  equation  (B-17a)  with  respect  to  x,  we 
obtain  an  expression  for  To 


(1  +  T0)6=  [b ExA(x')dx' . 
Jo 


(B-23) 


Now  multiply  equation  (B-21)  by  the  conjugate  of  the  aperture  field,  E*A, 
and  integrate  with  respect  to  x  using  equation  (B-22b)  to  find 
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(B-24) 
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Dividing  by  Jq  E*A(x)d.x  and  transposing  terms,  we  have 
I  §  +  if  jj <ix /»*•'««  (*,„  |.t  - 1'|)  £-a(x)£iA  (i') 
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Dividing  equation  (B-25)  by  equation  (B-23)  gives 
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The  input  admittance  at  z  =  0  is  defined  in  terms  of  the  reflection  coefficient 
[3,  p  14],  or 


i-r0 


(B-27) 


Yho  i  +  r0 

in  which  Yho  =  1/Zho  is  the  characteristic  admittance.  Using  equations 
(B-4)  and  (B-7)  with  n  =  0,  we  also  note  that  Znn  —  kznoJno/k2o  so  that 


Zhq  _  kzo 
Zhu  kzn 


(B-28) 


Substitute  equation  (B-22a)  for  Tn  in  equation  (B-26),  and  using  equations 
(B-27)  and  (B-28)  we  get 
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(B-29) 


which  is  a  variational  expression  for  the  input  admittance.  If  we  assume  a 
constant  aperture  field,  ExA{x)  =  Ex A,  then  the  first  term  on  the  right  side 
of  equation  (B-29)  vanishes  and  so  an  approximate  expression  for  the  input 
admittance  in  this  case  (i.e.,  Hio-excitation)  is 


^  “  it  [ dx  lo  dx'H°]  {kz0  |x  "  •  (B'30) 

In  terms  of  the  dimensionless  variable  x  =  kzox  —  kx  sin  9,  we  have 
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This  result  is  a  Green's  function  type  integral  in  which  the  integrand  de¬ 
pends  only  on  the  absolute  value  of  the  difference  parameter,  t)  —  \x~  X'\> 
so  is  an  even  function  of  77.  In  general,  we  are  interested  in  an  integral  of  the 
form,  I  =  Jo  Jo  f(x~  X')dxdx'  in  a  square  parameter  space  0  <  x  <  b  and 
0  <  x'  <  b.  Define  the  difference  parameter  asrj  =  x~x!  and  let  £  =  xj)  then 
as  shown  in  figure  B-3,  the  parameter  space  transforms  to  0  <i<b  and 
— £  <  77  <  b  —  £.  We  integrate  a  strip  over  77  in  region  1  and  region  2  where 
77  assumes  positive  or  negative  values,  respectively.  Referring  to  figure  B-3, 
the  double  integral  becomes 

=  /  Hv)dri  f  d^+  (  f  (77)^77  f  d£=  [  {f  (77)  +  /  (-77)}  (b  -  77)  dr] .  (B-32) 

JO  JO  J-b  J- 7?  JO 


Thus  if  f(?7)  is  an  even  function  of  r/,  then  equation  (B-32)  reduces  to 

1  =  2 b  [  f  {rj)dr]  —  2  [  r)f  {rj)dr) .  (B-33) 

Jo  Jo 

Now  with  b  =  kz0b,  let  7]  =  kz o{x  -  xf)  so  that  f(r/)  =  H ^  (kz o  \x  -  xt |),  is 
an  even  function  of  rj.  Then  with  equation  (B-33)  with  £o  =  1207 r  Q  the  free 


Figure  B-3. 
Transformation  of 
parameter  space  for 
Green's  function  type 
integrals. 
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space  impedance,  equation  (B-31)  can  be  written 
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For  cylindrical  Bessel  functions,  Z,,( Q,  we  know  that  f  C,Zv(QdQ  =  (,Z„+\  (£) 
and  so 
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Then  using  equation  (B-35)  to  simplify  the  second  term  in  equation  (B-34), 
we  find 


2k 

Yho 


kzo 


(kz0b)  +  lim 

77^0 


v  kz0b  Jj  ■ 


(B-36a) 


Where  we  use 


lim 

77 


kzob  J 


2  j 

irk-ob 


(B-36b) 


The  normalized  input  admittance,  ym,  for  the  parallel  plate  guide  radiating 
into  a  half-space  is  thus  from  equation  (B-36a)  with  equation  (B-36b) 


2/in 


H^(V)dV  -  H\ 


(2) 


(■ kzob )  +  j 


2 

Trkz()b 


(B-3  7) 


Expanding  the  Hankel  functions  into  Bessel  and  Neumann  functions, 
Hi  (()  =  J„  (£)  -  jN„  (C),  equation  (B-37)  is  written  in  terms  of  the  input 
conductance  and  susceptance,  normalized  to  the  iT-type  mode  admittance 


</in  —  v  —  f  Jo(v)dv  ~~  J\  ( kzob )  (B-38a) 

1  ho  Jo 
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(B-38b) 


B  rk*°b  2 

bin  =  vr-  =  -  /  N0(rj)di]  +  iVj  (/e20&)  +  -7- 
1  HO  Jo  7T  kz 

It  is  important  to  note  that  the  expression  (B-38b)  in  T/ie  Waveguide  Handbook 
includes  a  sign  error  in  the  first  term  [3,  p  184,  equation  (2a)].  Now  when 
the  aperture  is  not  the  full  waveguide  width,  d  <  b,  then  the  variational 
expression  can  still  be  used  to  obtain  equivalent  circuit  parameters.  The 
rationale  for  this  approach  is  that  with  constant  aperture  field,  the  equiv¬ 
alent  circuit  depends  only  on  the  external  field.  This  can  be  seen  from  the 
variational  expression  (B-29)  where  the  first  term  that  depends  on  the  inter¬ 
nal  fields  vanishes  when  EXj 4  is  constant.  Then  for  a  narrow  aperture,  the 
equivalent  circuit  parameters  are  given  by  equation  (B-38)  with  b  replaced 
by  d. 

The  normalized  circuit  parameters  in  equation  (B-38)  represent  the  "exact" 
solution  as  obtained  by  the  variational  method.  A  closed  form  approxima¬ 
tion  can  be  obtained  with  the  small  argument  formulas  for  the  Bessel  func¬ 
tions.  From  the  series  representations,  we  have 

Jo(*0  =  1  +  0  (z2)  ,  and  N0 (z)  =  ~  In  +  O  (z2)  (B-39a) 

Ji (*)  =  \  +  O  (z2)  ,  and  Ni(z)  =  ~  (j)  +  \ ln  (i)  +  “  (ln§  “  |)  +  0  (*2)  (B-39b) 


in  which  g  =  ec  =  1.781...,  with  the  irrational  number  c  =  0.577...  being 
Euler's  constant.  Thus,  when  kz0d  <  1,  equation  (B-38a)  becomes  (with  b 
replaced  by  d) 


9  in 


kzod 

~T~ 


kzod 

~2~ 


(B-40) 


This  small  argument  approximation  for  the  equivalent  conductance  is  the 
same  as  the  static  result  provided  by  Marcuvitz  [3,  p  184,  equation  (lb)]. 
The  equivalent  susceptance  from  equation  (B-38b)  becomes  (with  b  replaced 
by  d) 


2  fkz od  ( gri\ 

b'm  ~  /  In  (  —  ) 

7T  Jo  V  2  ) 


7 r 


kzod\  kzod 
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7T 


ln '  ^  —  2 


(B-41) 


Integrating  and  simplifying,  we  get 


kzpd  j  /  2e^/e  \ 
\gkz0d  J 


7 r 


(B-42) 


in  which  the  irrational  number  e  =  2.71828...  is  the  base  of  the  natural  log¬ 
arithm.  On  the  other  hand,  Marcuvitz  gives  the  static  solution  for  the  sus¬ 
ceptance  as  [3,  p  184,  equation  (2b)] 


kzod  . 

- In 

7 r 


(B-43) 
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which  cannot  be  obtained  from  the  truncated  series  representation  of  the 
Neumann  functions.  A  comparison  of  the  variational  result  for  the  suscep¬ 
tibility  equation  (B-38b)  to  the  static  result  equation  (B-43)  and  the  small 
argument  approximation  (B-42)  is  provided  in  figure  B-4.  The  small  argu¬ 
ment  approximation  equation  (B-42)  is  closer  to  the  variational  result  for 
kzod  <  0.5.  For  larger  kzod,  the  static  result  is  a  better  approximation  than 
equation  (B-42)  although  the  difference  is  negligible  for  k-od  <  1,  our  case 
of  interest.  Then  equation  (B-43)  is  used  when  a  closed  form  approximation 
is  required;  otherwise,  the  "exact"  susceptance  is  calculated  numerically 
from  equation  (B-38b). 


Figure  B-4.  Normalized 
susceptance  for  a 
parallel  plate  fed 
aperture  radiating  into  a 
half-space.  -g 
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Appendix  C.  The  Transverse  Resonance  Method 


The  electromagnetic  (EM)  fields  within  a  hollow  rectangular  waveguide  are 
represented  as  an  infinite  number  of  modes  that  are  propagating  or  attenu¬ 
ating  along  the  guide.  The  EM  field  of  each  mode  can  be  factored  into  trans¬ 
verse  and  longitudinal  field  components  with  respect  to  the  direction  of 
propagation  allowing  a  transmission  line  description  for  each  mode.  Even 
with  discontinuities,  a  waveguide  can  then  be  completely  described  by  an 
infinite  number  of  transmission  lines  [3,  p  6].  Following  Marcuvitz,  we  first 
summarize  the  conventional  (i.e.,  transverse  to  y)  mode  decomposition  for 
the  rectangular  waveguide  shown  in  figure  C-l  with  boundary  surface  S  in 
the  cross-sectional  plane.  Next,  we  consider  modes  transverse  to  z  with  the 
iT-type  mode  decomposition  derived  in  appendix  A.  We  place  the  origin 
as  shown  in  figure  C-l  so  that  a  transverse  transmission  line  analysis  in  the 
cross-sectional  plane  corresponds  to  propagation  along  the  positive  2-axis. 
We  describe  the  transverse  resonance  method  by  applying  the  technique  to 
the  closed  waveguide.  We  then  apply  the  method  to  a  localized  waveguide 
discontinuity  such  as  a  zero-thickness  aperture  in  the  waveguide  narrow 
wall. 


C.l  Rectangular  Waveguide  Modes 

We  only  consider  transverse  electric  (TE)  or  H -modes,  characterized  by 
modal  fields  with  Ey  =  0  while  Hy  must  satisfy  a  wave  equation  for  prop¬ 
agating  or  attenuating  waves  along  the  y-direction  with  longitudinal  prop¬ 
agation  constant  7  =  j(S.  When  completely  filled  with  a  uniform  dielec¬ 
tric  material  having  permittivity,  e  =  eo€r,  and  permeability,  y,  =  /uo/nr, 
the  modal  structure  is  the  same  as  for  the  air-filled  waveguide  but  with  a 
different  guide  wavelength  and  characteristic  impedance.  Referring  to  fig¬ 
ure  C-l,  take  the  .T2-plane  as  the  transverse  plane  and  define  the  transverse 


Figure  C-l.  Hollow 
rectangular  waveguide 
(with  perfect  condctor 
walls)  filled  with  a 
uniform  material. 
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(to  y)  Del  operator,  S7ty  =  +  J^z.  Then,  from  the  Maxwell  field  equa¬ 

tions,  the  EM  fields  must  satisfy 

=  M*  !  +  ^t!jJty  .  (H ty  x  y)  ,  and  (C-la) 

=  juj£  * +  "  •  (y x  E'«)  (c-u>) 

in  which  I  =  xx  +  yy  +  zz  is  the  unit  dyadic.  The  longitudinal  fields  are 
juyHy  =  Vty -(y  x  Eiy)  ,  and  (C-2a) 

juieEy  =  V(y  •  (Htj/  x  y)  (C-2b) 

This  representation  is  equivalent  to  the  Maxwell  equations  but  allows  sep¬ 
aration  of  the  transverse  and  longitudinal  field  variation  [3,  p  4].  In  the 
spirit  of  separation  of  variables,  expand  the  transverse  (to  y)  fields  in  vec¬ 
tor  modal  functions  that  are  constructed  to  have  the  correct  x  and  z  spatial 
variation  for  each  mode  i, 


Ety(x,y,z)  =  Vi(y)ei(x,z) ,  and  (C-3a) 

Htj,(x,  y,  z)  =  Ii(y)hi(x,  z)  (C-3b) 

in  which  i  denotes  a  double  index  mn.  The  vector  functions  define  the 
cross-sectional  form  of  the  mode  fields  while  the  voltage  and  current  func¬ 
tions  are  the  root  mean  square  (rms)  amplitudes  of  the  transverse  EM  fields 
at  any  point  y  along  the  propagation  direction.  Without  discontinuities, 
substituting  equation  (C-3)  into  equation  (C-l)  results  in  an  infinite  set  of 
transmission  line  equations 


and 


(C-4a) 


dlmn ( y ) 
dy 


j  hymn  Ymn  Vmn  (v) 


(C-4b) 


which  define  the  longitudinal  variation  of  the  mode  amplitudes.  Each  mnth- 
mode  corresponds  to  a  transmission  line  with  Zirm  =  1  /Ymn  the  character¬ 
istic  (wave)  impedance  and  jkymn  the  propagation  constant  along  y.  When 
we  use  equation  (C-3)  in  equation  (C-l)  and  take  account  of  equation  (C-4), 
the  transverse  (to  y)  mode  functions  are 


kymn  E rn n 


T  ,  YhiZll 

k2 


x  y)  , 


(C-5a) 


h,„„  =  [I  +  .  (y  x  e„)  .  (C-5b) 

Kymn*mv  L  & 

The  transverse  EM  fields  are  then  obtained  from  equation  (C-3)  in  terms  of 
the  transmission  line  voltage  and  current  functions.  With  equation  (C-3)  in 
equation  (C-2),  the  "longitudinal"  field  components  are 


juyHy  =  V7n„(y)Vty  ■  { yxe„m)  ,  and  (C-6a) 
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jusEy  —  Imn{y)^ ty  '  (hron  y)  • 


(C-6b) 


With  hmn  =  y  x  emn  and  emn  =  y  x  Vtyi>mn,  the  classical  TE-modal 
decomposition  is  obtained  in  terms  of  vector  modes  derivable  from  a 
scalar  mode  function.  The  mode  function  satisfies  the  scalar  wave  equation, 
(Vfy  +  k2mn)  i pmn(x,  z)  =  0  subject  to  the  boundary  condition  =  o,  on 

S  for  unit  outward  normal  vector  v.  Suitable  vector  modal  functions,  cho¬ 
sen  to  be  orthonormal  over  the  waveguide  cross  section. 


®mn*  Gpqdxdz  —  Snpdrnq  — 


for  n  =  p 
otherwise 


and  m,  =  q 


(C-6c) 


for  all  m  =1,  2,...  and  n  =  0,  1,  2,...,  allow  the  cross-sectional  dependence 
to  be  eliminated  from  equation  (C-l).  The  H-mode  functions,  emn,  normal¬ 
ized  over  the  waveguide  cross  section  are  then  derivable  from  the  scalar 
function 


‘Ipmn  0,  Z) 

with 


_  N/CmCn-^mn 1 


cos  (2ft)  COS  (: 
_  /  m2b  i  n2d 
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rrnrz  \ 
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M  V  a  1  ~b~ 

1  m  =  0  f  1  n  =  0 

2  m  ^  0  and^“\2n^0 


and  n  =  0, 1,2, . 


(C-7) 


in  which  Amn  is  the  mnth-mode  amplitude  normalization  factor  and  the 
time  dependence,  c]u:t ,  is  suppressed.  The  wavenumber  corresponding  to 
the  angular  frequency,  oj  is  written  in  terms  of  longitudinal  and  transverse 
components,  k2  =  ft2  +  k2mn.  Then  the  mode  function  satisfies  the  scalar 
wave  (Helmholtz)  equation  with  "cutoff"  or  transverse  wave  number 


=  *2  -  p2  =  (=) 


2  /  mr  \ ' 

+ 


T) 


(C-8) 


for  the  mnth-mode.  When  we  use  equation  (C-7)  in  equation  (C-5),  we  find 
the  vector  mode  functions  are 


/  x  -  ^  ,  VCmCnArm  fm n  .  / mnz\  „  m r  .  f  nnx\  A 

emn{x,  Z)  =  y  x  Vtytpmn  = - i== .  =  S - sin  -  x - —  sin  — —  I  z  1 

/m26  +  n2a(a  \  a  J  b  \  b  J  ) 


hm„(x,  z)=y  xev 


vCmCn^rrm  f  n7T  .  /mrx\  m7T  .  /m7TZ\  „ 

sm  — —  x  -| - sm  -  z 

V  b  J  a  \  a 


+  1  6 


and 

(C-9a) 

(C-9b) 


The  H -field  can  be  written  in  terms  of  the  (transverse)  mode  function  equa¬ 
tion  (C-7)  as 


tt  _  jPZo  t  t  \x  /  \  -  -  lCnn(y)  T-72  /  /  \ 

Hjnn  —  ~ :  Imn{y)R-mn\XiZ)  +  y  ;  V 

jup  jup  J 


(C-10) 
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and  xj>mn  satisfies  the  scalar  wave  equation,  so  the  EM  fields  are 

H)71n  =  ImT1(y)hmn{x,  z)  +  y-  7  z)  ,  and  (C-lla) 

W//  JUJfl 

E,)m {x.  y,  2)  =  %)  (C-llb) 

where  the  characteristic  impedance  is  the  ratio  of  the  voltage  to  current 
at  any  terminal  plane.  We  also  know  that  the  transverse  field  components 
propagate  down  the  guide  as  plane  waves.  Then  for  propagation  along  y 
we  must  have 


E  E 

y  x  =  z§nmn.  so  Zq  =  — —  =  — —  .  (012) 

il  x  jTZ  2 

Then  from  equation  (C-ll),  we  see  that  Zq  =  u)[i/(3  is  the  waveguide  charac¬ 
teristic  (wave)  impedance  for  //-modes.  The  EM  fields  are  thus  expressed 
in  terms  of  completely  transverse  vector  mode  functions  and  longitudinal 
modal  "transmission  line"  voltage  and  current  amplitudes.  The  voltage 
and  current  mode  amplitudes  defined  by  equation  (C-4)  satisfy  the  one¬ 
dimensional  wave  equation  so  are  either  propagating  or  attenuating  waves 
along  y.  For  propagating  waves,  an  impedance  description  using  stand¬ 
ing  waves  or  a  scattering  description  using  traveling  waves  can  be  used  to 
represent  the  voltage  and  current  mode  amplitudes  leading  to  a  complete 
description  of  the  modal  fields  [3,  p  9-13]. 

For  the  special  case  when  m  =  1  and  n  =  0  (i.e.,  the  Hio-mode),  the  funda¬ 
mental  (or  principal)  mode  EM  fields  from  equation  (C-ll)  simplify  to 

E10  =  -Ai()Vio(y))J~^sm  .  and  (C-13a) 

Hi° =  Ws  {iSVlo(!')  cos  (?)  * + ho(v)  sin  (?)  *}  e~il>"  ■ 

(C-13b) 

We  assume  that  only  the  fundamental  mode  propagates  so  that  from 
equation  (C-8)  the  transverse  wavenumber  is  A'no  =  7 r/a.  For  average  or 
rms  field  quantities,  power  flow  is  given  by  the  Poynting  vector,  Pjo  = 
Re  (E10  x  H£0)  .  The  cross  product  from  equation  (C-13)  is 

E,„  x  H;o  -  |^10|2  (1)  (H)y  -  |V„(,)|»  (§)’  i-n  (2)  cos  (=)  z}  . 

(C-14) 

The  mode  amplitude  is  determined  in  terms  of  the  power  flowing  at  some 
point  y'  within  the  guide.  The  input  power  Pjn(j/),  at  some  reference  plane 
y'  is  Pio(3’/,  y',  z')  integrated  over  the  waveguide  cross  section.  Since  equa¬ 
tion  (C-14)  is  independent  of  x,  we  have 

Pin(y')  S  b  r  P  10(?/,  z')d,z'  =  |^110|2  Re  (Vw(y')lM))  y  .  (C-15) 

J  0 
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The  fundamental  mode  amplitude  is  evidently  then  the  normalization  re¬ 
quired  to  have  |Pin(y,)|  =  Re(Vio(y,)Iio(y/))  at  some  reference  plane  y  =  y' , 
so  Aio  =  1  is  a  convenient  normalization.  The  transmission  line  character¬ 
ization  of  the  mode  is  not  unique,  and  different  normalization  constants 
could  be  used  [3,  p  8].  However,  absolute  impedance  comparisons  are  typ¬ 
ically  unnecessary  so  the  mode  description  can  be  left  in  terms  of  a  single 
amplitude  coefficient  determined  from  an  appropriate  normalization. 

We  now  summarize  waveguide  modes  that  are  transverse  to  the  z-direction, 
as  developed  in  appendix  A.  The  vector  EM  fields  E mp  and  Hmp  are  ex¬ 
pressed  in  terms  of  the  mode  function, 

ipp(x,  y)  =  Ap  cos  e~jl3y  (C-16) 

with  amplitude  normalization  factor  Ap.  This  mode  function  satisfies  the 
scalar  wave  equation 


(VL  +  ktP)  ^p(z,  y)  =  o  (C-17) 

with  transverse  (to  z)  wave  number  for  mode  p,  k2p  =  fc2  —  k2m,  and  Vte  = 
Now  the  wavenumber  is  decomposed  into  components  with 
respect  to  the  z-direction,  or 

k2  =  k%  +  k2zm  =  /92+(y)2  +  kL  ■  (C-18) 

The  propagation  constant  in  the  transverse  direction,  kzmr  is  unknown  but 
can  be  found  via  the  transverse  resonance  technique  by  defining  voltage, 
Vp(z)  and  current,  Ip(z)  mode  amplitudes  for  the  z-directed  transmission 
line.  These  functions  satisfy  transmission  line  equations  for  propagation  in 
the  transverse  direction,  or 

=  jkZmZHpIp(z)  ,  and  (C-19a) 

_dl^z)  =  jkgmYHpVp(Z)  (C-19b) 

CLZ 

where  the  characteristic  impedance  for  the  pth-mode  is  from  appendix  A 

z«v  =  WZtf  ■  <c-2°) 

The  vector  mode  functions  for  the  closed  waveguide  are  iT-type  modes  (see 
app  A) 

hP(x,  y)  =  -Vtzipp(x,  y)  ,  and  (C-21a) 

ep(z,  y)  =  h p(x,  y)  x  z  .  (C-21b) 
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Compared  to  equation  (C-2)  for  standard  H-modes  propagating  along  y, 
equation  (C-21)  yields  hybrid  modes  in  this  direction.  The  vector  EM  fields 
are 

E mp  —  Vp(z)ep(x,  y)  ,  and  (C-22a) 


H 


mp  -  Ip(z)hp(x,y)  +  z =  Ip(z)hp(x,y)  +  z -r^-il’P(x,y)Vp(z) 

JLOfl  JUJjl 


(C-22b) 


where  the  amplitudes  Vp{z)  and  Ip(z )  can  be  found  for  the  pth-mode  in 
terms  of  the  voltage  and  current  amplitude  at  a  single  point  [3,  p  9]. 


C.2  Uniform  Waveguide 

The  closed  waveguide  in  the  transverse  direction  can  be  considered  as  the 
junction  of  two  shorted  transmission  lines  as  shown  in  figure  C-2.  For  now, 
let  both  lines  have  the  same  characteristic  impedance  Zo  =  Zhp,  and  prop¬ 
agation  constant  yz  =  jkzm.  The  transverse  resonance  condition  requires 
that  at  resonance  the  input  impedance  (or  admittance)  "looking"  in  each 
direction  must  sum  to  zero  at  any  point  0  <  zt  <  a  [4,  p  167].  That  is 


Zm{z')  +  Zm {z')  =  0;  when  0  <  z'  <  a,  or  (C-23a) 

Ym  (z1)  +  Yin  (. z ')  =  0,  when  0  <  z'  <  a  (C-23b) 

where  the  arrows  denote  the  impedance  (or  admittance)  in  the  decreasing 
or  increasing  2-direction,  respectively.  At  a  point,  z'  located  a  distance  l 
from  the  load  impedance,  Zi(l),  the  input  inpedance  for  purely  imaginary 
propagation  constant  reduces  to  the  well-known  expression 


Z\n  (Y)  =  Z  (O  -  Zq 


ZliQ  +  jZp  tan  ( kztJ ) 
Z0  +  jZi(l)t&\\(kzml) 


(C-24) 


where  the  load  impedance  and  distance  to  the  load  on  the  left  or  the  right 
of  z'  are  to  be  used,  for  the  input  impedance  looking  to  the  left  or  right, 
respectively.  For  z'  =  0,  we  have  /  =  0  so  with  Z/(0)  =  0,  we  see  from 
equation  (C-24)  that  the  impedance  looking  to  the  left  vanishes.  Looking 
to  the  right,  the  termination  is  also  a  short  circuit  or  Z/(«)  =  0,  with  l  = 
a  -  z'  =  a.  Then  using  equation  (C-24)  with  Z0  =  ZHp  in  equation  (C-23a), 
we  find  the  resonance  condition  is  simply 


Zin(0)  =  jZHptan(kzma)  =  0  (C-25) 


which  requires  that  kzm  =  (mir/a).  Thus,  the  transverse  resonance  method 
applied  to  a  uniform,  closed  waveguide  (equivalent  to  a  transmission  line 
with  shorted  terminals  at  2  =  0  and  z  =  a)  yields  kzm  =  (mn/a)  as 
expected. 


48 


Figure  C-2.  Waveguide 
transverse  transmission 
line  description  as  a 
junction  of  two  shorted 
lines. 
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For  the  closed  rectangular  waveguide,  the  wavenumber  magnitude  from 
equation  (C-18)  is  then 

*=^T^=^+(^)2+(™)2.  (C-26) 

So  with  an  H- type  mode  decomposition,  the  transverse  resonance  condi¬ 
tion  results  in  the  Hmn-mode  waveguide  fields,  as  expected.  The  H- type 
hybrid  mode  fields  may  also  be  of  interest,  and  we  show  that  the  funda¬ 
mental  IT-type  mode  reduces  to  the  Hio-mode.  The  magnetic  vector  mode 
function  from  equation  (C-21a)  is 


hP0,  y )  =  -ApVtz  (■ 


COS 


p  KX 


*-jPy 


-*{(£)*  + 


pKX 


y  >e 


-jPy 


Then  using  this  in  equation  (C-21b),  we  find 


ep(x ,  y)  —  hp(x,  y)  x  z  =  Ap  j j/3  cos 


P'KX 


Pit 

x - —  sm 

b 


(C-27) 

(C-28) 


which  has  both  x-  and  y-components.  Compared  to  the  standard  TE-TM 
mode  decomposition,  these  functions  represent  hybrid  modes  having  both 
Ey  and  Hy  fields  in  the  longitudinal  (y-)  direction.  In  the  special  case  of 
p  =  0,  equation  (C-28)  becomes 

eo(z,  y)  =  A()j(3e~lliy-x. .  (C-29) 


For  a  shorted  transmission  line,  the  voltage  must  vanish  at  each  end  so  a 
Fourier  sine  series  could  represent  the  transverse  voltage 


0°  /  X 

Vp(z)  =  Y.  Brnp  sin  (  j  ■  (C-30) 

The  sum  is  over  a  possibly  infinite  number  of  modes,  each  with  amplitude 
Bmp,  a  complex  constant.  From  the  transmission  line  equation  (C-19a),  we 


have 

I  (2)  =  1 

dz  j  Zj{pkzm 


k‘2  -  (32  dVp(z)  1  _  .fc2-/32dyp(^) 

kzmuy  dz  jkzm  J  k2m up  dz 


(C-31) 
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so  that  the  pth-mode  current  is 

A»('4)  —3  ....  2^ 


w/<  “  \k'lma. 

r  rn  =  1  x 


m  7T 


Bmp  COS 


(=r)- 


(C-32) 


The  EM  fields  can  now  be  found  and  for  the  case  p  =  0  the  El-field  is 

E,„0  =  j/3e-^  ^  Bmo  sin  (  —  1  x  .  (C-33) 

m=l  \  ®  / 

The  normalization  constant  v40  is  absorbed  in  the  mode  amplitude  Bmn. 
The  iT-field  when  p  =  0  is 


H 


mO 


=  -  {  J:  £  (— )  B-.OCOS  (=i)  y  +  ^  £  ft*,  sin  ( z}  T-"*  • 

H'  ,^1  V  “  /  Vo/  w/*  V  o  )  j 

Then  for  the  typical  case  of  interest  m  =  1,  the  waveguide  fields  are 


(C-34) 


Eio  =  jfie  0sin  (  —  ]  x  .  and 

a 


TTZ 


Hio  =  —B\n  <  —  ( — 


a  \up 


cos 


TTZ\ 


jp2  .  ( KZ 


a  / 


y  +  - —  sin 

up 


z  >  e 


-jPv 


(C-35a) 


(C-35b) 


Eio  x  H\0  = 


These  H- type  fields  differ  from  the  H10-mode  fields  equation  (C-13)  by  a 
complex  constant.  The  amplitude,  I >\\>  is  determined  by  requiring  that  the 
average  power,  Pio  =  Re(Eio  x  (H*0),  integrated  across  the  waveguide  is 
the  input  power,  Pm{y').  The  cross  product  from  equation  (C-35)  is 

^  ^  |5io|2sin2  y  +  ^  |Bio|2sin  (^)  cos  z  (C-36) 


< up 


up 


which  is  purely  real  and  independent  of  x,  and  with  kto  =  P  from  equation 
(C-26)  we  have 


Pin(y')  =  b  fp  10dz  =(—)~  \B10\2y  (C-37) 

Jo  y  j  j 


The  normalization  constant  is  thus  Bi0  =  y  ^  so  from  equation  (C-35) 
the  EM  fields  are 


(C-38a) 

(C-38b) 


In  comparison  to  the  Hio-mode  fields  in  the  previous  section,  we  find  that 


^ioVjo(?/)  _  — j 
Moho{y')  -j 


(C-39) 


or  the  longitudinal  voltage  and  current  amplitudes  at  some  point  y'  are  re¬ 
lated  by  the  Hio-mode  impedance  Zq.  Thus,  the  fundamental  H- type  mode 
degenerates  to  the  Hio-mode. 
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C.3  Non-uniform  Waveguide 

The  transverse  resonance  method  can  also  be  applied  to  the  more  general 
case  of  a  rectangular  waveguide  filled  with  different  but  uniform  materials 
as  shown  in  figure  C-3.  The  closed  waveguide  in  the  transverse  direction  is 
now  a  junction  of  shorted  transmission  lines  having  different  propagation 
constants  and  characteristic  impedance.  As  for  the  uniform  guide,  H- type 
modes  are  used  for  the  transverse  mode  decomposition  in  each  region.  The 
characteristic  impedance  for  the  transmission  lines  along  2  to  the  left  and 
right  are  ZxHp  and  Z2Hp,  with  propagation  constant  71  zm  =  jhzm  and  7 2zm  = 
jkizmi  respectively.  We  outline  the  transverse  resonance  method  applied  to 
the  non-uniform  waveguide  but  omit  the  details  as  our  interest  here  is  in 
an  air-filled  waveguide. 

The  transverse  resonance  technique  that  uses  modes  transverse  to  the  z- 
axis  is  used  to  determine  the  transverse  wavenumber.  To  the  left  at  a  dis¬ 
tance  l  =  a0  away,  the  termination  is  a  short  circuit,  or  Z/(0)  =  0,  so  that  the 
input  impedance  looking  to  the  left  at  z'  =  ao  is 


Zin  (®o)  — 


£K0)  +  Zjfp  tanh  (jkizma0) 
ZlHp  +  Zi{  0)  tanh  (jklzma0) 


=  jZlHp  tan  (fci2ma0)  .  (C-40) 


The  transverse  line  to  the  right  is  also  terminated  in  a  short  circuit  at  a 
distance  l  =  a  —  ao  away,  or  Z/(a)  =  0.  Then  the  input  impedance  looking 
to  the  right  at  z'  =  ao  is 


7.  t  \  _  72  (a)  +  Zjjp  tanh  (jk2z 

a°  Hp  Z]ip  +  Zt  (a)  tanh  (jk2z 


(o  -  a0)) 


(a  -  ao)) 


=  j Zfjpt&n  (k2z 


(a  -  ao)) . 


(C-41) 


When  the  wavenumbers  are  purely  real,  the  resonance  condition  at  z'  =  ao 
requires 


jZffp  tan  (kizmao)  +  jZ%p  tan  (/c2zm  (a  -  ao))  =  0  .  (C-42) 


Figure  C-3.  A  closed 
waveguide  transverse 
description  as  a  junction 
of  different 
transmission  lines. 
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This  transcendental  equation  is  typically  written  in  terms  of  dimensionless 
variables  since  the  wavenumber  in  one  region  can  be  written  in  terms  of 
the  other.  One  equation  for  the  two  unknown  propagation  constants  is  ob¬ 
tained  by  subtracting  the  wavenumber  in  each  region.  Applying  the  trans¬ 
verse  resonance  condition  equation  (C-23a)  at  the  material  interface,  z'  =  0,0 
provides  another  equation.  Simultaneous  solution  then  determines  the  al¬ 
lowed  values  of  the  transverse  wavenumber  for  two  possible  cases,  k2 zm 
purely  real  and  k‘2Zm  purely  imaginary,  which  are  treated  separately.  The 
corresponding  waveguide  mth-mode  longitudinal  phase  constant,  fi,  then 
determines  the  phase  per  unit  length  in  the  propagation  direction  which 
completes  the  solution. 

For  the  nonuniform  line  created  by  geometric  changes  in  the  transverse  di¬ 
rection,  the  approach  is  slightly  modified.  The  discontinuity  is  lumped  at  a 
terminal  plane,  z  =  ao  +  h,  and  represented  by  a  discontinuity  reactance, 
jXo  as  shown  in  figure  C-4.  The  input  impedance  looking  to  the  right  is 
then  a  composite  impedance  formed  by  the  parallel  combination  of  the  dis¬ 
continuity  impedance  and  the  input  impedance  of  the  transmission  line  to 
the  right  terminated  in  load  impedance,  Z\.  The  input  impedance,  looking 
to  the  right  at  this  location,  is 


5  ln  ,  ,  x  _  V2  Z/(°)  +  ZHp.i  tan  {k'2zm  (o  -  0-0  -  /])) 

in(°°  °  HpZjIp  +  Z,(a)j  tan  (k2zm  (a  -  a0  -  /,))  ' 


(C-43) 


Then  at  z'  =  ao  the  transmission  line  to  the  left  is  terminated  with  the 
combined  impedance 


Z,(a  o  +  h)  =  ^^(«°  +  /i)  (C-44) 

jXd  +  Zlu  (ao  +  h) 

which  can  be  readily  calculated,  given  the  discontinuity  reactance.  The 
transverse  resonance  condition  applied  at  z'  =  oo  then  proceeds  as  before 
with  the  input  impedance  to  the  left  unchanged  but  to  the  right  we  have 


2  ,  .  _  ,71  Mao  +  ^i)  +  Zjjp, j  tan(kj;m/] ) 

P Z}ip  +  Z,(a0  +  l\)j  tan(kl2rn/i) 


(C-45) 


The  resonant  condition  then  yields  the  system  of  equations  to  solve  for  the 
transverse  wavenumber  in  each  region.  The  details  are  omitted,  but  the 
approach  as  outlined  again  demonstrates  the  versatility  of  the  transverse 
resonance  method. 


Figure  04.  Waveguide 
transverse  description 
with  a  discontinuity 
represented  by  a 
reactance. 
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C.4  Waveguide  Discontinuities 


The  transverse  resonance  method  can  also  be  applied  to  the  general  case  of 
a  rectangular  waveguide  having  a  localized  (i.e.,  zero  thickness)  disconti¬ 
nuity  in  the  transverse  direction.  The  discontinuity  could  be  an  aperture  in 
the  waveguide  narrow  wall  at  2  =  a  or  any  other  zero-thickness  structure 
that  can  be  represented  by  an  equivalent  circuit.  In  many  cases,  the  discon¬ 
tinuity  is  characterized  by  its  admittance.  Yd,  that  terminates  a  transmis¬ 
sion  line  with  characteristic  admittance  Yq  =  Yhp  for  each  of  the  p  possible 
modes.  The  transverse  resonance  condition  equation  (C-23b)  is  applied  to 
the  input  admittance  at  some  point  ao  as  shown  in  figure  C-5.  The  transmis¬ 
sion  lines  could  have  different  parameters  but  for  now,  we  let  Y^p  —  Y^p, 
and  71  zm  =  72 zm  =  jkzm ,  corresponding  to  a  uniform  waveguide  with  a 
localized  discontinuity. 

To  the  left,  the  termination  at  distance  l  =  ao  away  is  a  short  circuit,  or 
Yi  (0)  =  00,  so  that  the  input  admittance  looking  to  the  left  at  z'  =  a0  is 


^in  (ao)  —  Yjfp 


kj(0)  +  YHp  tanh  (jkzma0) 
Yhv  +  Yi(  0)  tanh  (jkzma0) 


jYfJp  cot  (fczm®o)  •  (C-46) 


The  transverse  line  to  the  right  is  terminated  at  distance  l  =  a  —  ao  away 
in  the  discontinuity  admittance,  or  Yj(a )  =  Yd-  Then  the  input  admittance 
looking  to  the  right  at  z'  =  ao  is 


^  /  '1  _  v  —£  -  Yjipj  tan  (fc  zm  (a  -  a0)) 

m  0  HP  YHp  +  YDj  tan  (kzm  (a -a0))  ’ 


(C-47) 


The  resonance  condition  equation  (C-23b)  applied  at  z'  =  ao  requires 


Yp  +  YHj)j  tan  ( k  zm  (a  -  a0)) 


...  n  \  .  -.t  -u  '  -  nPJ - r  "u;; 

-jYHp  cot  [kzmao)  +  YHp  - — — 1  - — - — - r- 

Yf{p  "1"  YdJ  tan  ( kzm  (a  ao)) 


=  0. 


(C-48) 


We  can  calculate  the  transverse  voltage  and  current  functions  for  mode  p 
by  finding  the  voltage  reflection  coefficient  at  2  =  ao-  To  the  right  for  mode 
P 

Vp(z)  =  A  sin  [kzm(z  —  a)]  for  ao  <  2  <  a  (C-49) 
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Figure  C-5.  Waveguide 
transverse  description 
with  a  discontinuity 
represented  by  a 
termination  admittance. 
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and  from  the  transmission  line  equation  (C-19a)  the  current  modes  are 

Ip(z)  =  jAYHp  cos  [kZ)„(z  -  a)]  for  a0  <  z  <  a  .  (C-50) 

The  transmission  line  on  the  left  is  terminated  at  z  =  a, a  with  the  input 
admittance  looking  to  the  right  equation  (C-47).  The  reflection  coefficient 
to  the  right,  for  each  mode  p  is  then 


fp(«o)  =  ^ . -  (C-51) 

The  voltage  and  current  to  the  left  of  ao  include  reflected  waves  from  2  = 
o0/  or 


Vp  (. 2 )  =  B  +  r p  (ao)  ejk^{z~ao) 


for  0  <2  <  a, o  (C-52) 


since  the  reference  is  at  z  =  ao,  where  rp(oo)  is  known  for  each  mode.  The 
pth-mode  current  is 

IP(z)  =  B YHp  {e-jkzmiz~ao)  -  rp(o0)e^'--”'(2-°o)}  for  0  <  ^  <  a0  .  (C-53) 

The  boundary  condition  at  z  =  ao  requires  continuity  of  the  current  and 
voltage  or 


A  sin  [kzw(a,o  —  a)]  —  B  {1  +  r(o0)}  =  0  ,  and  (C-54) 

jA  —  B  {1  —  T(o0)}  =  0  .  (C-55) 

Nontrivial  solutions  of  this  system  of  equations  are  possible  only  if  the  de¬ 
terminate  (det)  of  coefficients  vanishes,  or  det  =  0.  This  condition  has  al¬ 
ready  been  satisfied  by  the  transverse  resonance  condition,  and  the  condi¬ 
tion  for  resonance  can  be  derived  by  setting  det  =  0  in  the  above  system.  The 
final  result  for  the  voltage  and  current  modes  depends  on  a  single  constant 
(the  mode  amplitude)  so  let  5=1,  then  from  equation  (C-55) 

A  =  —  ^(ao)  =  j  (T(ao)  -  1)  (C-56) 
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The  voltage  and  current  amplitude  for  each  mode  p  are  then  determined; 
thus,  the  structure  of  the  waveguide  fields  can  be  evaluated. 

We  are  interested  in  the  case  when  ao  =  a,  with  fundamental  H- type  mode 
propagation  in  the  transverse  (z-)  direction  or  p  =  0.  So  for  an  Hio-excited 
waveguide,  we  use  m  =  1  and  have  a  reflection  coefficient  at  the  plane 
containing  the  discontinuity  (e.g.,  an  aperture)  given  by 


Yho  -  fin  _  YEl  ~  ft  __  Ye 0  ~  VD 
Ym  +  Ym  Ym  +  Yg  Yho  +  vd 


(C-57) 


in  which  yn  is  the  normalized  discontinuity  admittance.  The  resonant  con¬ 
dition  equation  (C-48)  becomes 


—j  cot(kzia)  +  =  0  ,  or  (C-58a) 

Yho 

hi  =  ~  cot-1  (-jyo)  (C-58b) 

a 

which  would  require  numerical  solution.  The  voltage  and  current  ampli¬ 
tudes  to  the  left  are  the  transverse  mode  amplitudes  throughout  the  wave¬ 
guide,  or 

V0(z)  =  j  (T  -  1)  sin  [kz i  (z  -  a)]  ,  and  (C-59a) 

Io(z)  =  (1  -  T)  Yho  cos  [kzl(z  -  a)]  =  (1  -  f )  —  cos  [kzl(z  -  a)]  (C-59b) 

(jJfJL 

with  T  =  To  (a)  as  given  by  equation  (C-57)  and  kz i  =  \J k7  —  j32,  given  by 
equation  (C-58b).  The  fundamental  mode  fields  throughout  the  waveguide 
with  localized  discontinuity  at  z  =  a  are  then 

Eio  =  A0(l  -  r)/3sin  [kzl{z  -  a)]  e~j/3y±  ,  and  (C-60a) 


Hio  =  A0  (1  -  T)  Y0  {jkz i  cos  [kz\(z  -  a)]  y  +  (5 sin  [kz\(z  -  a)]  z}  e  jl3y  (C-60b) 

in  which  To  =  P/wp  is  the  Hio-mode  characteristic  admittance  and  Ao  is 
the  mode  amplitude.  Of  course,  these  EM  fields  are  perturbed  near  the  dis¬ 
continuity,  but  for  a  "small"  discontinuity  or  at  sufficiently  large  distances, 
the  perturbation  would  be  negligible.  This  represents  the  complete  solu¬ 
tion  as  obtained  by  the  transverse  resonance  method  for  this  special  case  of 
interest. 
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Appendix  D.  Numerical  Routines 


The  numerical  calculations  presented  were  accomplished  with  MATLAB™. 
The  routines  as  implemented  in  the  form  of  "m-files"  are  included  for  com¬ 
pleteness.  Minor  modifications  of  these  files  may  be  required  to  produce  all 
the  results  shown.  The  subroutine  for  the  FEM  calculations  is  included  at 
the  end  of  the  main  routine. 


o  /o  /o  /o  /o  /o  /o  / o  /o  /o  /o  /o  /o  /o  /o  fo  fo  /o  /o  fo  /o  /o  /o  /o  / o  / o  /o  fo  /o  / o  /o  /o  / o  /o  /o  /o  yo  /o  yo 


%  NOTAPER.m  (CONSTANT  WIDTH  SLOT) 


'0  /0  /O  /0  /O  /O  /O  /O  /O  /O  /O  /O  /0  /O  /0  /O  /O  /0  /0  /0  / 0  /O  /O  /O  / 0  70  /o  /o  /0  70  70  70  70  70  70  70  70  70  70 


%FOR  PROPAGATION  IN  WG  PERTURBED  BY  NARROW- WALL  SLOT 

%COMPUTES  COMPLEX  PROPAGATION  CONSTANT  VS.  DISTANCE. 

%TE10-MODE  PROPAGATION  IN  THE  INPUT  AND  OUTPUT  WAVEG¬ 
UIDE 


'O  /0  /0  /0  /0  /0  /O  /0  /O  /0  /0  /0  /0  /O  /O  /0  /O  7 0  /0  /0  /o  /0  /o  70  70  70  70  fo  70  70  70  70  70  70  70  70  70  70  70 


%  METHOD  OF  SOLUTION: 


%COMPLEX  PROPAGATION  CONSTANT  BY  PERTURBED  TRANSVERSE 
WAVE 


%FOR  NARROW  SLOTS  PERTURB  ABOUT  SHORT  CIRCUIT  (see  ECE 
239  5/2000). 


0  /e  /o  /o  /o  /o  /o  /o  70  /o  /o  70  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /0  /0  70  70  70  70  70  70  70  70  70  70  70 


%FEM  FOR  INPUT  REFLECTION  (TRANS)  COEFF,: 

%PROVIDES  ANTENNA  EFFICIENCY  AS  1  -  R2  -  T2  FOR  REALIZED 
GAIN 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  REMARKS: 

%  WAVEGUIDE  CROSS-SECTION:  a  =  WIDE  DIM.  ALONG  X, 
%HALF-HEIGHT  WG:  HEIGHT:  b  =  NARROW  WALL  ALONG  Y  (E-aperture) 
%NARROW-WALL  SLOT  WIDTH:  d  ;;  b/10;  PROPAGATION  ALONG  Z 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  PLOT  OUTPUTS: 

%  1)  SLOT:  ALPHA  &  CURRENT  VS  NORMALIZED  DISTANCE 
%  2)  PHASE:  BETA/BETA10  VS  NORMALIZED  DISTANCE 
%  3)  FIELD:  FARFIELD  PATTERN  (H0  and  //G) 

%  4)  RADIATION:  REALIZED  GAIN  (H-  and  C-plane) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  2)  PERTURBATION  SOLUTION  FOR  NARROW  SLOTS  ONLY 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  WOC/ 27  JULY  2000 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


tic; 


%MM=input('Input  number  of  FEM  cells  per  guide  wavelength: '); 


MM=100;  delta=le-10; 
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% . Free  space  parameters.. 


e=exp(l);  g=1.781; 
mu=4e-7*pi;  c=3e8; 

% . Waveguide  parameters . 

a=2;  b=a/ 2; 

%  SET  WAVENUMBER  FOR  OPERATING  FREQUENCY 

k0=4.25/  a;  lam=2*pi/k0; 

omega=kO*c; 

lgl0=lam./sqrt(l-(lam/2/ a).2); 
betal0=sqrt(k0.2-(pi/ a)2); 

Z10=120*pi*lgl0.  /  lam; 

%  Standard  sphereical  coordinates  theta  from  z  and  phi  from  x 

%Current  along  y-axis  with  endfire  angle  theta=phi=pi/2 

Nth=200;  Nph=100; 

thmax=60*pi/ 180;  phmax=pi; 

dph=phmax/Nph;  dth=thmax/Nth; 

th=delta:dth:thmax-delta; 

ph=delta:dph:phmax-delta; 

% . Slot  parameters . 


'Input  uniform  slot  maximum  width 


dO=input('normalized  to  wg  height : '); 


dO=dO*b; 

%%%%%%%SLOT  LENGTH  NORMALIZED  TO  WG  WIDTH 
'Input  length  of  slotted  section  in  wavelengths 
L=input(['Total  slot  length  normalized  to  wavelength  : ']); 
L=L*lam; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%-/0707b 


'0/..0L 


M=MM*ceil  (L/lglO); 

N=M+l;dy=L/(N); 

y=delta-L/2:dy:L/2; 

%%%STEP  1:  Uniform  slot  width  vs  axial  position 
d=dO*y./y;  dmax  =  max(d); 

%%%STEP  2:  Narrow  slot  width  perturbation  formulation  at  each  d(y) 
G=pi.*d/a/2; 

B=(b*log(csc(pi.*d/2/b))+d.*log(a*e/ g./ d))/a; 

A=G.2+B.2;  R=G.  /  A;  X=B./A; 


kxO=pi/ a  +  B./A/a  +  j.*G./A/a; 


beta0=sqrt(k02-kx0.2); 

%% TEST  RESULT  TO  ENSURE  NEGATIVE  ATTENUATION  CONSTANT. 


for  mm=l:N; 


if  imag(beta0(mm)));0;  betaO(mm)=betaO(mm)-2*i*imag(betaO(mm)); 
end;  end; 

alpha=-imag(betaO); 

%%%STEP  3:  Call  FEM  code  for  solution  with  TE10  at  each  end 
betaO(l)=betalO;  betaO(N)=betalO; 

%%%%%H-type  hybrid  mode  impedance  (see  ECE297  3/ 99) 
ZH=omega*mu./real(betaO); 


V  /O  /0  /O  /0  /O  /0  /0  /O  /O  /O  /O  /C  /O  / O  /0  /O  /0  /0  /O  /O  /O  /O  /O  /0  /0  /O  / 0  /O  /0  /0  /O  /O  /O  /O  /O  /O  /O  /O 


[phi]=onedfem(beta0,M,dy,Z10,ZH); 


'0  /O  /O  /O  /0  /O  /O  /O  /0  /O  /O  /O  /0  /O  /O  /O  /0  /O  / 0  /O  /0  /O  /O  /O  /O  /O  /O  /O  /O  /O  /O  /O  /O  /O  /O  0  /O  /O  /O 


%%%STEP  4:  CALCULATE  REFLECTION,  TRANSMISSION  &  ATTENU¬ 
ATION 


RR=((phi(l))-l); 

TT =((phi(N))*exp(i*betaO(N)*L)) ; 
prad=l  -  abs(RR)2  -  abs(TT)2;, 
eff  =  ceil(prad*100); 


%%%STEP  5:  CALC  AMPLITUDE  DISTRIBUTION  FROM  ALPHA 
for  mm=l:N; 

temp(mm)=alpha(mm);  aO=alpha(mm); 
F2(mm)=a0*exp(-2*dy*sum(temp)); 
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end;  M0=sqrt(F2); 


%%% STEP  6:  CALCULATE  FAR-FIELD  FOR  UNIFORM  ILLUMINATION 
F(y)  =  1 

for  n=l:Nth;  for  m=l:Nph; 

ARG=kO*sin(th(n))*sin(ph(m))-real(betaO); 

I2=abs(dy*sum(M0.*exp(i*y.*ARG)))2; 

%  Hth2(n,m)=cos(th(n))2*sin(ph(m))2*I2; 

%  Hph2(n,m)=cos(ph(m))2*I2; 

H2(n,m)=cos(th(n))2*sin(ph(m))2*I2... 

+  cos(ph(m))2*I2; 
end; 

%%%%%%%%%%%% Factor  of  four  included  for  integration  limits 

phiint(n)=dph*sum(H2(n,:))/  pi; 

end; 

denom=dth*sum(phiint.*sin(th)); 

%%%%%STEP  7:  CALCULATE  DIRECTIVE  GAIN  (prad  FOR  REALIZED 
GAIN) 

Gp=10*logl0(2*L/lam); 

['Expected  gain  (2L/lamda) is  num2str(Gp,3)] 

%  %  %  %  PATTERN  FROM  NUMERICAL  INTEGRATION  OF  CURRENT%%%% 
D  =  prad*H2/denom;  clear  H2; 


G  =  10*logl0(D);  clear  D; 


maxG=max(G(:,Nph/  2)); 


for  n=l:Nth; 

if  G(n,Nph/ 2)==maxG;  bwO=n;  end;  end; 

['Integrated  peak  gain  (in  dBi)  is  num2str(maxG,4)] 


V  /o  /o  /o  /o  /o  70  /o  /o  /o  / 0  fo  /o  /o  70  /o  /o  70  70  /o  70  70  70  /o  fo  /o  /o  /o  /o  fo  /o  /o  /o  /o  70  fo  /o  /o  /o 


th=th*180/ pi;  ph=ph*180/ pi; 
amax  =  max(alpha  /  kO) ; 
bmin  =  min(real(betaO)/betalO); 

%%%% . plot  output . 

figure(l); 

subplot(2,l,l),  plot(y /L,alpha*lam); 
axis([-.5  .5  0.8*amax,  1.2*amax]); 

textf-.S,  1 . 1 5*amax,  ['Slot  radiation  efficiency  = ',  num2str(eff)/ 
text(-.25,1.05*amax,['  kOa  = ',  num2str(k0*a)/  with  a  =  ',num2str(a),'  m.']); 

text(-.2,0.9* amax, ['Total  slot  length  = ',  num2str(L/lam,3)/  wavelengths.']); 

title('Slotted  waveguide  alpha  (in  Nepers/ wavelength).'); 

ylabel('Normalized  loss  (alpha*lamda)') 

subplot(2,l,2),  plot(y /  L,M0  /  max(MO)); 

text(-.25,l/ 2, ['Maximum  slot  width  = ',  num2str(dmax,3), '  m. ']); 
text(.3,l/2,['L  =  ',  num2str(L/ a)/a ']); 
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title('Current  amplitude  versus  normalized  distance.'); 


xlabel('Normalized  distance  (y/L)'); 
ylabel('Normalized  current7); 
figure(4); 

plot(th,G(:,Nph/2)/k-'); 

xlabel('Elevation  angle,  theta  (degrees)'); 

ylabel('Realized  gain  (in  dBi)'); 

title('H-plane  gain  for  a  slotted  waveguide  antenna.'); 

text(25,l.l*maxG,['Numerical  realized  gain  =  num2str(maxG,4),'  dBi.']); 

text(35,maxG,['On  boresight  phi  =  ',  num2str(ph(Nph/2),3)/  degrees.']); 

figure(5); 

plot(ph,G(bwO,:)); 

xlabel(' Azimuthal  angle,  phi  (degrees)'); 

ylabel('Realized  gain  (in  dBi)'); 

title('E-plane  gain  for  a  slotted  waveguide  antenna.'); 

text(25,l.l*maxG,['Numerical  realized  gain  = ',  num2str(maxG,4),'  dBi.']); 

text(35,maxG,['At  main-beam-angle,  thetaO  = num2str(th(bw0),3),'  degrees.']); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


['Total  clock  time  is ',  num2str(toc/60,3),'  minutes.'], 
%ONEDWAVE  is  a  non-uniform  t-line  code  using  FEM  as  in  Jin 


%  INPUT  AXIAL  PROPAGATION  CONSTANT  AND  FEM  PARAMETERS 


%  4/ 00  function  [phi]=  onedfem(beta0,M,dz,Z10,ZH); 

%Uses  onedfem  to  solve  for  TE  wave  propagation  (Jin  p.54) 

%for  non-uniform  T-line  in  axial  direction  w/o  sources  so  f=0 
%Using  M  finite  element  cells  and  linear  interpolation  N=M+1; 

%Element  lengths  constant  mesh  l=ones(M,l)*dz; 

%BVP  parameters  for  non-uniform  T-line  propagation  alpha=-l./betaO./ ZH; 
beta=betaO./ZH; 

%B.C.  from  T-line  eqs  for  reflection  and  transmission  at  terminal  planes 
gammaL=-i/ZH(N);  qL=0; 
gammaO=-i  /  Z10;  q0=-2*i  /  Z10; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%Solves  a  1-D  BVP  using  FEM  for  M  elements  using  Jin  notation 

%Solution  (phi=V  or  I)  Subject  to  mixed  (or  Neuman)  B.C.  at 

%at  z=0  -alpha*dE/  dx+gamma*E=q  and  at  z=L  alpha*dE/  dx+gamma*E=q 

%Setup  interaction  matrix  size  NxN 

K(l)=alpha(l)/l(l)  +  beta(l)*l(l)/3; 

K(N)=alpha(M)  /  1(M)  +  beta(M)*l(M)  /  3; 
for  n=2:N-l; 

K(n)=  alpha(n-l)/l(n-l)  +  beta(n-l)*l(n-l)/3  +  ... 
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alpha(n)/l(n)  +  beta(n)*l(n)/3; 


end 

for  m=l:N-l; 

C(m)=-alpha(m)/l(m)  +  beta(m)*l(m)/6; 
end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


for  k=l:N;  b(k)=0;  end; 

%Enforce  Mixed  (or  Neuman)  B.C.  at  x=0 
K(l)=alpha(l)/l(l)  +  beta(l)*l(l)/3  +  gammaO; 
b(l)=qO; 

%Enforce  Mixed  (or  Neuman)  B.C.  at  x=L 
K(N)=alpha(M)/l(M)  +  beta(M)*l(M)/3  +  gammaL; 
b(N)=qL 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%Solve  system  for  V(z):  slash  uses  Gaussian  Elimination  for  k  =  2:N; 


K(k)  =  K(k)  -  ( C(k  -  1  f/K[k  -  1)); 


b(k)  =  b(k)  -(C(k-l)*b(k-l)/ K(k-l)); 


end  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


phi(N)=b(N)/K(N); 


for  p=l:N-l; 
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k=N-p;  phi(k)  =  (b(k)  -  (C(k)*phi(k+1)))/K(k); 
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